AD-A046  293 


GTE  SYLVANIA  INC  NEEDHAM  HEIGHTS  MASS  ELECTRONIC  SYS— ETC  F/G  9/3 
CPMFSK  AS  A MODULATION  ON  THE  VLF-HF  ATMOSPHERIC  CHANNEL. (U) 

OCT  77  T A SCHONHOFF  F30602-76-C-0410 


N 


UWf.LASSIFIFn 


RADC-TR-77-318 


References 


1.  R.  R.  Anderson  and  J.  Salz,  "Spectra  of  Digital  FM" , BSTJ , 

Vol.  44,  pp.  1165-1189,  July-August  1965. 

2.  W.  P.  Osborne  and  M.  B.  Luntz,  "Coherent  and  Noncoherent  De- 
tection of  CPFSK",  IEEE  Trans.  Comm. , COM-22,  #8,  August  1974, 
pp.  1023-1036. 

3.  T.  A.  Schonhoff,  "Symbor  Error  Probabilities  for  M-ary  Co- 
herent Continuous  Phase  Frequency  Shift  Keying",  ICC  1975 
Conf . Reed. , pp  34-5  to  34-8. 

4.  T.  A.  Schonhoff,  "Symbol  Error  Probabilities  for  M-ary  CPFSK: 
Coherent  and  noncoherent  Detection",  IEEE  Trans.  Comm., 

COM-24,  #6,  June  1976. 

5.  J.  T.  Gamble,  "Linear  Detection  of  Coherent  FSK  (CFSK)  and 
MSK  Signals  in  Low  Frequency  Atmospheric  Noise",  ICC  1975 
Conf.  Reed. , pp.  34-14  to  34-18. 

6.  "World  Distributions  and  Characteristics  of  Atmospheric  Radio 
Noise",  CCIR  Report  322,  1964,  Figure  27. 

7.  G.  D.  Forney,  Jr.,  "The  Viterbi  Algorithm",  Proc.  IEEE,  Vol. 
61,  No.  3,  March  1973,  pp.  268-278. 

8.  A.  D.  Whalen,  Detection  of  Signals  in  Noise,  Academic  Press: 
New  York,  1971. 

9.  J.  N.  Craig,  Design  of  Lossy  Filters,  M.I.T.  Press,  Cambridge, 
Mass.  1970. 

10.  A.  J.  Viterbi,  "Error  Bounds  for  Convolutional  Codes  and  an 
Asymptotically  Optimum  Decoding  Algorithm",  IEEE  Trans.  Int. 
Thy. , IT-13,  April  1967,  pp.  260-269. 

11.  T.  A.  Schonhoff,  "Bandwidth  vs.  Performance  Considerations 
for  CPFSK",  NTC  1975  Conf.  Reed. , pp.  38-1  to  38-5. 

12.  H.  M.  Hall,  Jr.,  "A  New  Model  for  'Impulsive'  Phenomena: 
Application  to  Atmospheric-Noise  Communication  Channels," 

Ph.D.  dissertation,  Stanford  University,  August  1966. 

13.  A.  Giordano  et  al.,  "Interference  Mitigation  Techniques 
Program",  interim  report  under  contract  N00014-76-C-0468 
sponsored  by  ARPA. 

14.  J.  W.  Modestino  and  A.  H.  Levesque,  private  correspondence. 

15.  O.  Ibukun,  "Structural  Aspects  of  Atmospheric  Radio  Noise  in 
the  Tropics,"  Proc.  IEEE,  Vol.  54,  No.  3,  March  1966,  pp. 
361-367. 


87 


16.  J.  W.  Modestino,  "Locally  Optimum  Receiver  Structures  for 
Known  Signals  in  Non-Gaussian  Narrowband  Noise",  13th  An- 
nual Allerton  Conf.  on  Ckt.  and  System  Theory,  Oct.  1975. 

17.  I.  S.  Gradshteyn  and  I.  M.  Ryzhik,  Tables  of  Integrals, 

Series,  and  Products,  Academic  Press:  New  York,  1965. 

18.  W.  Q.  Crichlow  et  al.,  "Determination  of  the  Amplitude-Prob- 
ability Distribution  of  Atmospheric  Radio  Noise  From  Statis- 
tical Moments",  Journal  of  Research,  National  Bureau  of 
Standards  - D:  Radio  Propagation,  Vol.  64D,  1960,  pp.  49-56. 

19.  K.  E.  Wilson,  "Analysis  of  the  Crichlow  Graphical  Model  of 
Atmospheric  Radio  Noise  at  Very  Low  Frequencies",  M.  S. 

Thesis,  Air  Force  Institute  of  Technology,  Nov.  1974. 

20.  A.  D.  Spaulding,  C.  J.  Roubique,  and  W.  A.  Crichlow,  "Con- 
version of  the  Amplitude-Probability  Distribution  Function 
for  Atmospheric  Radio  Noise  from  One  Bandwidth  to  Another," 
Journal  of  Research,  National  Bureau  of  Standards  - D Radio 
Propagation,  Vol.  66D,  No.  6,  Nov. -Dec.  1962. 

21.  J.  T.  Gamble,  "An  Analysis  of  Linear  and  Nonlinear  Coherent 
Detection  in  Atmospheric  Noise  At  Very  Low  Frequency",  Ph.D 
dissertation,  Syracuse  University,  1974. 

22.  A.  Shohara,  "VLF/LF  Modem  Performance  in  Atmospheric  Noise- 
Volume  1:  Theoretical  Analysis",  Interim  Final  Report, 

January  1976,  DCA  100-75-C-0013. 

23.  C.  D.  Beach  and  D.  C.  George,  "Error  Performance  of  VLF  and 
LF  Receiving  Systems  with  Nonlinear  Atmospheric  Noise  Reduc- 
tion", RADC-TR-70-190 , Westinghouse  Electric  Corp. , Sept. 

1970.  (AD875991) 

24.  M.  G.  Pelchat,  R.  C.  Davis,  and  M.  B.  Luntz,  "Coherent  Demodu- 
lation of  Continuous  Phase  Binary  FSK  Signals",  1972  Int. 
Telemetering  Conf.,  pp.  181-190. 

25.  0.  Y.  Antonov,  "Optimum  Detection  of  Signals  in  Non-Gaussian 
Noise",  Radio  Eng.  and  Electronic  Physics.  V.  12,  April  1967, 
pp.  541-548. 

26.  S.  S.  Rappaport  and  L.  Kurz,  "An  Optimal  Nonlinear  Detection 
for  Digital  Data  Transmission  Through  Non-Gaussian  Channels", 
IEEE  Trans.  Comm.,  COM-14,  No.  3,  June  1966,  pp.  266-274. 

27.  M.  Abramowitz  and  I.  A.  Stegun,  Handbook  of  Mathematical 
Functions,  National  Bureau  of  Standards,  AMS  55. 

28.  A.  Papoulis,  Probability,  Random  Variables  and  Stochastic 
Processes , McGraw-Hill,  New  York:  1965 , Section  7.2. 


88 


Appendix  A 

Selected  Theoretical  Results 

This  Appendix  contains  details  of  the  theoretical  effort. 
The  emphasis  here  is  the  general  determination  of  the  probability 
of  error  when  the  received  signal  has  been  processed  in  a bandpass 
nonlinearity  followed  by  a linear  receiver.  The  basic  approach 
that  has  been  chosen  is  to  model  the  output  of  the  linear  receiver 
as  a finite  sum  of  random  variables  (note:  they  need  not  be  inde- 
pendent nor  identically  distributed).  As  explained  in  Section 
5.1.2,  we  can  define  a variable  . (Eq.  78) 

L 

y = E a(y,  (Al) 

l=o  * * 

where 

y £ = C(t£)  cos  (C(t£)  - 6 ( t £ ) (A2) 

All  the  terms  are  defined  in  Section  5.1.2.  We  will  use  a type  A 
Gram-Charlier  Series  to  model  the  pdf  of  y. 

Gram-Charlier  Series 

Other  authors  have  used  the  Gram-Charlier  approach  to  at- 
tempt to  determine  the  performance  of  modulation  systems  in  im- 
pulsive noise  when  nonlinearities  are  being  used  [22,23].  Beach 
and  George  [23]  used  a severely  truncated  (three  or  four  terms) 
series  and  their  results  do  not  match  empirical  results  [21]  and 
Shohara  [22]  uses  an  approach  which  requires  that  all  the  samples 
be  independent  and  identically  distributed.  This  is  applicable 
for  PSK  modulation,  but  does  not  lend  itself  to  CPMFSK  demodula- 
tion. Shohara  recognized  this  and  indicated  in  an  Appendix  ways 
which  the  problem  might  be  solved. 

The  derivation  of  the  Gram-Charlier  series  is  given  in 
many  references  [8,  22,  23].  Thus,  just  an  overview  will  be  given 
here.  If  we  define  a new  variable 
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z = 


Y-Wi 


0 (A3) 

where  y = E(y}  and  a = Var  {y},  the  pdf  of  z can  be 

given  as  [8] 


p(z)  = 4>  ( z)  + I 

£=3  bi<t>{z)  H£(z) 


(A4) 


where 


1 z 

<f>(z)  = — — . e “2 


(A5) 


and 


,th 


H (z)is  the  £ Hermite  polynomial. 


H^(z)  can  either  be  determined  by  a recursion  formula  [8] 


H£(z)  = zH£-1(z)  " (H)H(  ,(z) 


(A6 ) 


where  H (z)  = 1 
o 

and 


(A7) 


H^(z)  = z,  or  it  can  be  shown  that 


|J.I+1  , ,,k  , n-2k 

H. (z)  - j i:.1)  Si  z 
k=o  (n-2k)!k!2k 


(A8) 


i £ i 

where  indicates  the  nearest  integer  less  than  or  equal  to 


To  determine  the  coefficients  b?  in  Eq.  A4 , we  can  multi- 
ply both  sides  of  A4  by  Hm(z)  and  integrate  over  the  real  line. 

(A9) 

00  oo  00  00 

/ H (z)p(z)dz  = / <M  z)  H ( z)  dz  + £ b f <f>  ( z ) H0  ( z)  H (z)dz 

-oo  1,1  .oo  m J,=  3 * _oo  £ m 
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As  long  as  m^l,  the  first  term  on  the  r.h.s.  of  A9  is,  from  [27, 
Section  22],  0.  Furthermore,  the  Hermite  polynomials  are  bior- 
thogonal  so  that 


1 


= mi  6 , 


Where  6n  is  the  Kronecker  delta, 

Jem 


6n  = 
£m 


1,  £=m 

0,  £/m 


Combining  A9  and  A10  we  get 


f H (z)p(z)dz=bm!  or 
— oo  in  m 


(A10) 


(All) 


(A12) 


b = -f-7—  f°°  H ( z ) p ( z ) dz 

mm!  m 

—00 


(A13) 


Since  Hm(z)  is  a polynomial  in  z,  knowledge  of  the  moments  of  z 
will  allow  the  determination  of  b , and  hence  from  A4 , p(z).  From 
A3,  we  know  that 


So  E {z  } = 


m 

/ TV>  \ 

1 

k) 

0 

11 

\ k / 

1 

m 

m 

y. 

a 

0 

1 n 

(k)  (-u  )k  E {ym-k} 


(A14) 


(A15) 


Hence,  if  we  can  determine  the  moments  of  y,  we  can  calculate  the 
moments  of  z.  The  next  section  discusses  techniques  to  determine 
the  moments  of  y. 
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From  Eq.  A4 , it  is  clear  that  a practical  determination 
of  p(z)  will  necessitate  a truncation  of  the  infinite  series.  Ni 
known  techniques  have  been  found  for  determining  the  truncation 
error,  and  this  is  a question  still  to  be  addressed. 

Before  concluding  this  section,  it  is  recalled  from  Sec 
tion  5.1.2  in  the  text  that  the  desired  calculation  is  to  deter 
mine  the  probability  that  y exceeds  zero.  This  implies,  from 
Eq.  A3,  that  we  desire 

Pr{z  > -^•1}  = f p(z)dz  (A16 

"yl 

a 

This  can  be  seen  to  be 


Pr{z  > — 1 } = / <f>(z)dz  + £ b /°°  4i(z)H^(z)dz  (A17 
-U  i £=3 

a a 


1 erfc  ( ]L I)  . (A18 

2 /2a 


From  [8] , we  know  that 


From  L 2 7 ] , 


o 


/27T 


-z 

2 


dz  = 


<P(z)Hz(z)  = (-1)  [ £]  (z) 

where  <f>^(z)  is  the  derivative  of  <Mz)  . 

/ (-1)  % U1  (z)  =(-^4>U~11(z)  = 


Therefore , 


(Al  9 


(-1) 


24-1 


$ (z)H  (z)  , (A 


so 


/.  <Mz)H#  (z)dz  = 


ZH.1 

a 


H 


/TrT  £-1 


(-Ui  ) 


y i 

2a2 


(A21) 


Hence 


“Vi 

Pr  { z > } 


1 

7 erfc  (-7=r— ) 

2 /?  a 


+ E b£H£ 
£=  3 ^ ^ 


-y- 


-yi 

2a2 


(A22) 


This  implies  that  it  is  unnecessary  to  use  the  Gram  Charlier 
series  to  determine  the  pdf  of  z,  but  it  can  be  used  directly  to 
determine  the  probability  that  z exceeds  a threshold. 


Determination  of  Moments  of  y 

From  Eq.  Al , it  can  be  shown  that 


(A23) 


k 

y = 


k ! 


m 


m. 


^ / _ \ w / 

mo=0  m^=0  mL=0  m0 ' oyo'  'alyl) 


nu 


'*(aLyL) 


with  the  added  condition  that 


L 

I m = k ( A24 ) 

£=0  * 

Eq.  A23,  which  is  a multinomial  expression  can  be  thought  of  as 
an  extension  of  the  binomial  theorem.  From  A2  3,  the  kth  moment 
of  y is 
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IP 


E {yK } = Z 


k 

Z 


m0=0  irij  =0 


k 

Z 


k! 


m^O  mo  ! nij  ! • • •mL' 


(A-25) 


m, 


E { y. 


in  r i 

• yT  L> 


with  the  condition  of  Eq.  A24.  If  knowledge  is  available  or  de- 
rivable about  the  joint  statistics  of  the  various  time  samples, 
then  Eq.  A2  5 can  be  used  to  determine  the  kt'1  moment  of  y.  On 
the  other  hand,  if  each  of  the  time  samples  is  independent,  then 

E{y>y^  = E{yom°}E{Ylmi}  ...  E{yLmL}  (A26) 


Eq.  A26  can  be  used  in  Eq.  A25  to  again  calculate  the  k*"*1  moment 
of  y,  however  in  the  case  of  independent  sample,  a recursive  tech- 
nique can  be  devised. 

From  Al 


yk  = (ll0  a*y*)k  = (a°y°+  aiyi+  * 


• + aLyL) 


( A27 ) 


Using  the  binomial  theorem  this  can  be  written  as 


»L  ("o)  laoyo,k'n°  ,aiy 


1 + a2y2  + • 


+ aLyL}  ° 


(A28 ) 


Using  the  binomial  theorem  on  the  second  term  of  A28  gives 


(A29) 


n*-o(»o)  ("0  '*lyl,B#  ni“*Y*  +-  + VlT1 
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and  repeated  application  gives 


k 

y = 


J.  W 


, . k-no  "o 

(a0y0)  z 

nl“0 


/no\ 

\nl/ 


(alyl) 


no-ni 


n 


L-2 


I 


L-2 


nL-2  nL-l 


nL-l=° 


^aL-lYL-l ^ 


L-l  (A30 ) 


n 


ULyI.1 


L-l , 


From  this,  taking  the  expected  value  of  both  sides  allows  a re- 
cursive  technique  for  calculating  E{y  }.  Clearly  either  tech- 
nique requires  the  determination  of  moments  of  the  individual 
samples  y^  , thus  the  next  section  will  address  the  determination 
of  these  momemts. 


Determination  Moments  of  y^ 


From  (A2) , 


y2,  = COS  ~ 6 (fc£) ) 


{ A31 ) 


where  from  the  text,  £(t^)  is  the  sampled  signal  plus  noise  en- 
velope which  has  been  processed  by  the  nonlinearity,  z,  (t^)  is 
the  sampled  signal  plus  noise  phase,  and  6(t^)  is  the  sampled 
filter  phase.  5(t^)  an<^  c(t^)  are  ran<^om  variables  and  6 ( "t ^ ) 
is  deterministic. 

Using  the  analytic  noise  models  of  Section  4,  the  joint  pdf 
of  £(t^)  and  C(t^)  can  be  derived  as  follows.  We  will  temporarily 
drop  the  sampled  representation.  Consider  figure  Al , which  shows 
a signal  vector  of  magnitude  A and  phase  ip  and  an  additive  noise 
vector  of  magnitude  n and  phase  <J> . We  are  then  interested  in 
the  resultant  signl  plus  noise  vector  of  magnitude  v and  phase  a. 
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I 

I 


SIGNAL  IS  Re{Ae  '*e  'wt} 

NOISE  IS  Re  i<f>e  iwt  J. 

4403-77E 


Figure  Al.  Envelope  and  Phase  of  Signal  Plus  Noise 
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It  can  be  seen  that 


2 ? i 

v = [ (A  cos  ip  + n cos  (p)  + (A  sin  1(1  + ri  sin  $ )‘1],  (A- 32) 


= [A^  + + 2An  cos  ( cf>  — ) ] 15 


and 


a = tan 


~1  A sini/;  + 


n sm 


A cosip  + n sin  <J> 


( A—  3 3) 


We  know  the  joint  pdf  of  n and  <p  , designated  as  p„  ,(  n , <P  ) . To 

H , 

find  the  conditional  joint  pdf  of  v and  a conditioned  on  A and  ip , 
we  need  the  Jacobian  of  the  transformations  (A32)  and  (A33)  [28]. 

With  simple  substitution,  the  Jacobian  is  n/v.  As  an  example, 
consider  the  truncated  Hall  model  with  0=3.  From  section  4,  the 
joint  pdf  of  n and  $ is 


PH , <J> ( n ' ^ ' = 


by^n 


u [n2+y2] 2 
0 


0 < n < e 


n > e 


(A-34) 


It  follows  that  the  conditional  joint  pdf  of  v and  a is 


pV , a (v»a I 


by  v [v2  + A^  - 2Av  cos  (a-ijj)]*8 
it  [v^  + A2  - 2Av  cos  (a-ij/)  + y 2]2 


(A-35) 


, 0 < n < E 


0,  n > E 
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The  signal  plus  noise  sample  passes  through  the  zero-memory  non- 
linearity which  is  modeled  in  Section  2 by  the  transformation  of 
variables . 


e = 


£ - a 


v,  v < th 


k,  v > t. 


(A- 36) 


where  k = 0 for  a bandpass  hole  punch  and  k=t^  for  the  bandpass 
clipper.  Using  (A-36),  the  joint  conditional  pdf  of  ?,  ? is 


P-'Z  U#C|  = wth  (?)  PVfA  (v,a|A,iH 


+ 6(?-k)  tf  pv  (v|A,t|;)dv  . pA  ( ? | A , ij; ) dv 

h 


( A-  3 7 ) 


where  wt^  (?)  is  a window  function  and  6 ( ?— k ) is  the  Direc-delta 
function  correspondingly  defined  as 


wth  (^> 


and 


5 (?-k) 


1-  C£th 

0,  ?>th 

1,  ?=k 

0,  elsewhere 


(A38) 


( A-39 ) 


p (v  A,^)  and  p (?|A,i|;)  are  the  conditional  marginal  densities 

V A 

for  the  envelope  and  phase  respectively.  Eq.  A-37  can  be  under- 
stood by  recognizing  the  first  term  as  that  part  of  the  input  pdf 
which  is  unaffected  by  the  nonlinearity  (input  envelope  is  less  than 
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the  threshold) , whereas  the  second  term  is  the  part  which  is  af- 
fected by  the  nonlinearity. 

It  is  important  to  note  that  the  joint  conditional  pdf  of 
£ and  t,  as  shown  in  Eq.  A-37  is  non-zero  on  a finite  range,  i.e. 
for  ^t^.  This  will  greatly  facilitate  anv  numerical  techniques 
which  must  be  used  in  further  processing. 

Now  we  make  one  further  change  of  variables 


y - £ cosU-6) 
x = C sin  U-<5)  , 


(A-40 ) 


the  joint  conditional  pdf  of  x and  y is  p (x,y|A,i|>) 

X , X 


(A-41) 


Py  y (*/y  |A,iW  P-  7 / y o 0 —1  v 

/ x2+y2  (/x2  + y2f  tan  (£)+6|a,i|») 


from  which  the  marginal  density  of  y is 


PY(y  I A , ) 


px  y (x,y |A,^)dx 


(A-42) 


The  mfc^  conditional  moment  of  y is  then 

E { ym  1 A , } = /th  ym  py  (y|A,U>)  dy  (A-43) 

_th 

In  all  the  foregoing,  we  have  suppressed  the  sampled  representa- 
tion. None  of  the  noise  models  known  to  the  author  allow  analytic 
determination  of  all  the  above  equations.  It  is  interesting  to 
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note  that  even  Gaussian  noise,  for  which  the  input  pdf  of  the 
envelope  is  a Rician,  does  not  allow  full  analytic  solution  to 
all  the  above  equations.  Hence,  numerical  approaches  must  be 
used  in  the  above  integrations.  This  in  turn  requires  a study 
of  the  resolution  required  for  the  numerical  integrations.  Time 
constraints  precluded  this  study  and  the  ensuing  numerical  deter- 
minations of  the  moments  of  individual  samples. 
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APPENDIX  B 


Statistical  Relationships  for  Hall  and  Truncated  Hall  Models 

This  appendix  contains  the  derivations  of  the  expressions 
presented  in  Table  5. 

Envelope  PDF  for  Truncated  Hall  Distributions;  9 = 2,  3 

The  envelope  pdf  for  the  truncated  Hall  distributions 
(0  = 2,  3)  are  defined  to  be 


PE  (n) 


0-1 

b (0-1) y n 

, (0  + l)/2 

(n  +y  j 


= 0 

where  b is  chosen  so  that 


0 < n < E 


n > e 


(B-l ) 


E 

/ p (n)dn  = l (b-2 ) 

0 t 

Using  the  transformation  n = ytanz,  0 < z < tan_1(E/y)  in 
Equation  (A-2)  leads  to 


0-1 


(B-3 ) 


where 


1+ (E/y) 
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Thus,  for  0=2,  the  envelope  pdf  is  given  by 


P (n)  = I D 

E D-l  (n2+Y2)3/2  ' 


= 0 


< n < yJd2-i 


, n > yJD  -i 


(B-4 ) 


and,  for  0=3, 


P (n)  = _°L  _2lLJ± 

E d2-i  (n2+y2) 2 


= o 


° < n < yJo2-i 
n > y Jo2-! 


(B-5 ) 


Moments  of  Envelope  pdf 


and 


The  moments  of  the  envelope  are  defined  by 


w i 

u-  = f n p ( n ) dp 


= r 19.7.1) -y9"^^1  dn  e = 4,  5 


2.  2 


0+1 


(n  +Y  ) 2 


(B-6 ) 
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D f (6-1)  Y n- 

0-1  J ^ ^ 9 + 1 

D - 1 0 (n  +Y  ) 2 


9-„Vi 


dn  0=2,3 


(B-7) 


The  transformation  used  to  obtain  ( B—  3 ) above  yields  expressions 
for  the  convergent  moments*. 

Thus  integrals  of  the  form 


7,  2'  2.  0+1  dn  ~ 0-j-l  f 

(n  +Y  ) -5-  Y J J 


j+1  2 

tan  z sec  2 
0+1 

sec  z 


fVr  /si 


j+1  0 — j — 2 

in  Z cos  z dz 


can  be  evaluated  for  specific  values  of  j and  0,  using  the  ap- 
propriate  limits  of  integration.  Substitution  of  D = Vi+e2/y2 
in  the  case  of  truncated  distributions  results  in  the  expressions 
given  in  Table  5.  For  example,  the  mean  of  the  envelope  pdf  for 
0 = 2 is  given  by 


"1  ' D?I  if 


E 2 


f -tv- = ^ 

0 (n  +Y  ) 2 


sin  z 


A*  { 


log  (D  + %/  D - 1 ) - 


■H 


*For  0 = 4,  moments  Uj,,  u2  converge;  for  0 = 5,  j = 1,  2,  3 

converge;  all  moments  of  the  truncated  distributions  exist. 
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Amplitude  Probability  Distributions:  APD 

The  amplitude  probability  distribution  provides  the 
probability  of  the  variable  exceeding  a level  R,  i.e.. 


-I 


Pr{n>R)  = / p(n)dn 


For  the  Hall  distributions  with  0=4,5 


Pr{ n>R>  = (0-l)y 


V 


/ 2 ? 9+1 
R (n  +Y  ) 2 


(B-10 ) 


(B-ll ) 


and  for  0=2,  3 


Pr{n>R)  = -L-  (0-1)y9" 

D - 1 


1 f ndn 

J 2 2 — 

R (n  +y  ) 2 


(B-12 ) 


2 2 

Using  the  substitution  z = n +y  / 


(B-13 ) 


(9'11  f 7T1*1  = <6'1)  /-if r - -5^r  - 

J (n  +Y  ) 2 j 2z  2 z 2 (n  +Y2)  2 
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Therefore  the  Hall  distributions  have  the  following  APD's: 


For  0=2 


Pr(n>R)  = Y 


(E2+y2)?  (R2+y 2}2  J 


D 

D-l 


For  0=3 


Pr{n>R)  = — — Y2 


For  0=4 

Pr { n>R)  = — 
(R 


-1 


(E2+Y2) 


2 2 
(R  +Y  ) 


(B-14 ) 


(B-15) 


(B-16) 


: 
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For  0 = 5 


Pr { n>R)  = 


(r2+y2) 


('*?) 


(B-17) 


Marginal  Density  of  Quadrature  Components 

The  amplitudes  of  the  quadrature  voltage  components  are 
related  to  their  envelope  v by  the  transformation 


n = n cos 
c 


n_  = n sin  4) 


(B-18 ) 


where  <(>  is  assumed  to  be  uniformly  distributed  and  independent 
of  n . 

Therefore,  for  0=4  and  5 where 


P(P,)  ( 2 ^ (e+1)/2 

(n  +y  J 


o < n < 


(B-19) 


the  transformation  (A-18)  gives 


oo 

iO  -/ 


(0-1) 

2tt 


2 22 
(nc  +n8  +Y  > 2 


(B -20) 


The  quadrature  components  n and  n are  identically  distributed 

c s 

and  p(ng)  = P(nc) • 


r 


1 


2 2 2 

Letting  a = nc  +Y  in  (B-20)  and  making  the  substitu- 
tion n = a tanz,  yeilds  integrals  of  the  form 
s 


/ 


ns 


2.  2 


0+1 


(ng  +a  ) 2 


1 f 6-1 

= ^0  J cos 


z dz 


for  which  recursion  formulas  exist.  For  0 = -5 


(B-21) 


, x 3y 

P'V 2 , 5/2 

4 (nc  +y  ) 


- 0°  < V 


(B-22 ) 


For  0=4 


P(nc} 
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2 2 
(nc2+Y2) 


(B-23) 


For  0=2  and  3,  the  marginal  density  is  obtained  by  substituting 


■ 


Appendix  C 

Statistical  Relationships  for  the  Generalized  Laplace  Variable 

The  probability  that  the  generalized  Laplace  variable,  s 
exceeds  the  threshold,  R,  is  given  by 


Pr { s>R}  = / p(s)  ds 

R 


(C-l) 


= 1 - 


R 

f 


s r K 


l-r<£> 


r+1  r-1 

(ko)  r (r) 2 


ds 


(C-2 ) 


Using  Item  6.561-12  in  [17]  to  evaluate  the  above  equation  yields 


Pr {s>R} 


= (k°)  r (r)  2r_1  K (*«) 


(C-3) 


Since  K 


/.L  \ (A] 

-r  \ko)  = Kr  ^ka J , 


equation  (45)  follows. 


If  s denotes  the  envelope  of  the  quadrature  components 
x and  y,  the  pdf  of  x (or  y)  can  be  found  using  the  trans- 
formations 


and 


x = s cos  <j> 
y = s sin  <p 


(C-4 ) 


where  <p  is  a uniformly  distributed  random  variable  in  the  interval 
(0,  2n) . 


L. 
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The  pdf  of  x can  then  be  found  from 


00 

p(x)*  = / p (x,y ) dy 

— oo 


J0  JL 

_o6  2tl 


(/x2+y2) 


r-1 


r+1 


(ko) 


T (r ) 2 


_ Ki_r(45ZL\dy 

r-1  \ ko  / 


(C-5) 


Noting  that  p ( — x)  — p (x)  and  using  item  6.596  — 3 in  [17]  yields 


p(x) 


r+*s  r-*s 

(ko)  itT  (r ) 2 
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Appendix  D 


Transformations  for  Probability  Graph  Paper 

The  transformation  of  the  abscissa  values  for  the  proba- 
bility graph  paper  as  used  here  disagrees  somewhat  with  the  trans- 
formations used  by  Crichlow  and  Wilson.  The  point  of  disagree- 
ment is  the  value  of  the  constant  multiplying  the  right  hand  side 
of  the  equation.  In  order  to  clarify  any  problems,  here  is  the 
derivation  of  the  abscissa  transformation,  when  the  ordinate 
transformation  is  known: 

y*  = 20  log1()y  (Dl) 

when  the  abscissa  transformation  is  assumed  to  be  of  the  general 
format : 


x'  = C log1Q  (-In  x)  (D2) 

and  assuming  that  the  Rayleigh  function  plots  as  a straight  line 
of  slope  -0.5.  The  Rayleigh  pdf  is 

2 

p(z)dz  = exp  (-  — - — ) dz  ( D 3 ) 

° 2 2a2 

Since  the  rms  value  of  z is  /2  a,  let  x = — , (04) 

/2  a 

2 

resulting  in  p(x)dx  = 2x  exp(-x  ) dx  (D5) 

The  APD  of  R is  APD  (R)  = /“  p(x)dx 

R 


= exp  (-R2) 

The  slope  of  the  Rayleigh  line  is  -0.5 


20  log  R 

C log (-In (APD(R) ) ) 


( D6 ) 


(D7) 
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20  log  R 
2C  log  R 


From  this  we  can  see  that 


C = -20 

and  Eq.  55  in  the  text  follows. 


( D8 ) 
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EVALUATION 


This  effort  investigated  the  quantitative  value  of  using  Constant  Phase 
Multi-Frequency  Shift  Keyinq  (CPMFSK)  modulation  and  multi-baud  demodulation  in 
impulsive,  Non-Gaussian  noise  typical  of  sub-HF  radio  channels.  Previous  work 
has  shown  that  CPMFSK  modem  techniques  are  much  more  efficient  than  more 
conventional  modems  in  Gaussian  noise,  and  also,  that  non-linear  impulsive 
noise  reduction  increases  in  effectiveness  as  the  detection-element's  length 
increases.  Thus,  the  combination  of  noise  reduction  and  CPMFSK  multi-baud 
detection  was  exoected  to  provide  enhanced  performance  in  real-world  atmos- 
pheric noise.  However,  the  results  of  this  study  indicate  that  the  relative 
advantaqe  of  CPMFSK  multi-baud-detection  (over  MSK)  and  non-linear  noise 
reduction  diminishes  with  increasing  noise  impulsiveness.  Thus,  the  attract- 
iveness of  CPMFSK  for  impulsive  VLF-HF  is  somewhat  reduced.  Another  aspect  of 
this  work,  an  investigation  of  noise  clipping  and  blanking  thresholds, 
reaffirmed  previous  assumptions  that  a hard-limiter  was  equal  or  superior  to 
any  other  bandpass  clipper  or  blanker  in  reducing  impulsive  noise;  however, 
an  equally  interesting,  and  potentially  valuable,  finding  was  that  a blanker 
is  significantly  more  effective  than  a clipper  when  the  threshold  is  constrained 
to  be  higher  than  the  mean  envelope  of  the  siqnal  plus  noise. 

This  effort  has  provided  definitive,  although  lower  than  originally 
conjectured,  performance  results  for  CPMFSK  signals  in  a realistic  sub-HF 
radio  environment.  Nevertheless,  the  findings,  in  predetection  receiver 
structures,  demodulation  efficiency,  and  signal  spectra,  definitely  contribute 
to^the  state  of  knowledge  in  radio  communication  technology. 

-tc 

^ JOHN  T.  GAMBLE 
Project  Engineer 
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1.  Introduction 


I 


The  VLF-HF  Atmospheric  Noise  Channel  is  characterized  by 
an  impulsive  like  noise  environment.  Much  of  the  present  digital 
communications  on  this  channel  employ  constant  envelope  modula- 
tions, and  most  of  these  modulations  can  be  considered  forms  of 
Frequency  Shift  Keying  (FSK) . Specific  types  of  FSK  modulations 
which  have  been  used  on  the  Atmospheric  Noise  Channel  are  Minimum 
Shift  Keying  (MSK)  and  Coherent  Frequency  Shift  Keying  (CFSK) . 
Both  of  these  modulations  can  be  described  as  specific  examples 
of  Continuous  Phase  Multi-Frequency  Shift  Keying  (CPMFSK) , i.e., 
the  modulation  uses  an  FSK  type  signal  wherein  the  phase  is  con- 
strained to  be  continuous  from  one  symbol  to  the  next.  This 
constraint  of  phase  continuity  can  have  two  beneficial  results 
on  the  qualities  of  the  signal: 

(1)  The  spectral  occupancy  of  the  signal  can  be  con- 
siderably smaller  than  traditional  modulations  [1]« 

(2)  The  signal  has  memory,  and  utilization  of  this 
memory  by  observing  the  received  waveform  for  more 
than  one  symbol  interval  can  result  in  significant- 
ly improved  performance  [2-4]  . 

Most  of  the  performance  improvements  determined  pre- 
viously have  assumed  an  additive  white  Gaussian  noise  (WGN ) 
channel.  It  is  interesting  to  note  that  for  the  WGN  channel, 
the  error-rate  performance  of  MSK  is  the  same  as  the  performance 
of  Phase  Shift  Keying  (PSK) . Yet,  reference  [5]  indicates  that 
on  the  VLF  Atmospheric  Noise  Channel  the  performance  of  MSK  is 
superior  to  that  of  PSK.  This  seems  to  be  due  to  the  fact  that 
MSK  demodulation  uses  essentially  2-bit  time  intervals  of  signal 
to  make  decisions,  whereas  PSK  (and  most  traditional  modulations) 
use  independent  symbol-to-symbol  decisions. 

In  [2-4]  it  is  shown  that  a more  general  CPMFSK  signal 
can  achieve  performance  results  which  exceed  MSK  in  the  WGN 
channel.  This  improvement  performance  is  due  largely  to  the  fact 
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that  more  memory  is  imposed  on  the  signal,  and  observation  inter- 
vals of  many  symbols  (n>^3)  are  utilized  to  make  decisions.  These 
facts  were  the  primary  motivations  for  the  present  study. 

This  report  documents  the  result  of  this  study  to  inves- 
tigate the  performance  of  CPMFSK  when  used  as  a modulation  on  the 
additive  Atmospheric  Noise  Channel.  It  was  the  intent  of  the 
study  to  determine  the  error-rate  performance  of  CPMFSK  by  both 
theoretical  estimation  and  Monte-Carlo  simulation.  Furthermore, 
the  effect  of  pre-detection  filtering  and  the  general  concern  of 
receiver  design  was  addressed.  Finally,  transmitted  signal  power 
spectral  densities  (p.s.d.)  were  derived  for  the  modulations  of 
interest . 

In  the  attempt  to  derive  theoretical  estimates  of  error- 
rate  performance,  it  was  found  necessary  to  derive  and  evaluate 
mathematically  tractable  Atmospheric  Noise  models.  These  models 
were  constrained  to  match  some  of  the  empirically  determined 
first-order  statistics  of  the  Atmospheric  Noise  Channel  [-o]. 

All  of  the  goals  of  the  study  were  achieved  except  for 
the  theoretical  estimation  of  error-rate  performance.  However, 
error-rate  performance  by  means  of  simulation  was  obtained.  This 
report  will  document  the  present  status  of  the  theoretical  effort 
but  time  and  budgetary  constraints  precluded  the  full  development 
and  results  of  this  phase  of  the  study. 

Section  2 below  outlines  the  mathematical  development  and 
the  definition  of  terms  which  should  greatly  augment  an  under- 
standing of  the  rest  of  the  report.  Section  3 considers  more 
realistic  receiver  implementations,  with  specific  emphasis  being 
given  toward  the  effects  of  various  types  oi  predetection  filters 
Section  3 also  considers  a receiver  implementation  based  upon  the 
Maximum  Likelihood  Sequence  Estimation  (MLSE^  algorithm,  or  as 
it  is  perhaps  more  commonly  known,  the  Viterbi  Algorithm  (VA) 

[7].  Section  4 present  the  noise  modeling  results  which  were 
accomplished  during  the  study.  Three  noise  models  were  used,  two 
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for  the  theoretical  effort,  and  one  for  the  simulation  effort. 
These  models  are  shown  to  adequately  match  first-order  statistics 
of  the  Atmospheric  Noise  Channel.  The  primary  results  of  the 
study  are  presented  in  Section  5.  A brief  review  of  the  theoret- 
ical effort  is  given,  however,  as  mentioned  above,  this  phase  of 
the  effort  was  not  completed,  so  the  bulk  of  the  theoretical 
effort  has  been  relegated  to  Appendix  A.  Section  5,  in  turn, 
consists  primarily  of  simulation  results  which  nonetheless  do 
indicate  that  CPMFSK  modulation  yields  a performance  improvement 
over  MSK , i.e.,  CPMFSK  can  achieve  error  rates  equivalent  to  MSK 
with  a smaller  signal-to-noise  ratio.  Section  6 presents  the 
power  spectral  densities  (p.s.d.)  of  the  CPMFSK  signals  of  inter- 
est. 
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2.  Mathematical  Background 
2.1  Signal  Representation 

th 

A received  CPMFSK  signal  during  the  1 symbol  interval 
can  be  written  as 


wct+ 

- 


d . 7rh  ( t-  (i-1 ) T ) 

1 s 


T 

s 


+ 


i-1 

I 

j=l 


irhdj 


+n  (t) 
(1) 


where  E is  the  symbol  energy,  T is  the  symbol  interval,  cj  is 
s s c 

the  carrier,  h is  the  deviation  ratio,  and  d^  is  the  digital 
information.  The  term  1 j;"''  nhdj  is  the  phase  at  the  beginning 
of  the  ifc^  interval  whi^h^can  be  assumed  zero  for  i=l.  n(t)  is 
additive  impulsive  noise,  which  can  be  written  in  a bandpass 
representation  [8]. 


n(t)  = n (t)  cos  (o  t -n  (t)  sin  w t 
c c s c 


(2) 


where  nc(t)  and  ng(t)  are  lowpass  impulsive  noise  processes. 

For  M-ary  modulation,  d^  = + 1,  + 3,  . . . , + (M-l)  and  any  two 
adjacent  frequencies  in  the  signal  set  are  separated  by  h/T  Hz. 


For  the  commonly  u^ed  MSK,  h=0.5  and  for  CFSK,  h=1.0. 
This  results  in  orthogonal  signaling  for  both  cases.  If  h is 
chosen  other  than  0.5  or  1.0  and  the  signal  is  observed  for  n 
symbol  intervals,  better  performance  than  orthogonal  signaling 
ran  result  and  performance  improvements  over  traditionally  used 
modulations  is  possible.  The  procedure  to  improve  performance 
is  to  observe  the  received  signal  for  n symbol  intervals,  and 
then  make  a decision  on  one  of  the  n symbols.  In  the  case  of 
coherent  detection,  the  decision  is  made  on  the  first  symbol. 


In  [2-4],  the  optimum  receiver  for  high  SNR  detection 
in  a Gaussian  noise  environment  was  derived.  It  is  useful  to 
use  the  shorthand  notation  developed  in  [2]  for  the  received 
signal.  


r (t) 


A (t,  d^ , Dk)  + n (t) 


(3) 
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where  is  the  first  symbol,  and  is  the  (n-1)  tuple 
Dk  = (d2, . . . , dnl . 

The  vector  can  of  course  take  on  Mn  1 = m different  possibil- 
ities. With  these  definitions,  the  optimum  receiver  for  high 
SNR  can  be  depicted  as  shown  in  Figure  1.  That  is,  the  receiver 
correlates  the  received  signal  with  each  of  the  Mn  possible 
n-symbol  waveforms  that  could  have  been  transmitted.  After  n 
symbol  intervals,  a decision  is  made  on  the  first  symbol. 

It  will  be  useful  in  the  sequel  to  use  a baseband 
equivalent  representation  of  the  received  signal  rather  than  the 
bandpass  representations  of  Eqs.  1 and  3.  Correspondingly 

j id  t 

d^ , D^)  = Re  { s ( t , d1#  Dk)e  } (4) 

Using  Eqs.  1,  3 and  4 we  can  see  that  the  baseband  equivalent 

th 

signal  waveform  during  the  i symbol  interval  is 


s (t,  d^ , D^)  = exp 


d.irh(t- (i-l)T  ) 


f rdiTih(t- 
\3L T 


i-1  -i 


+ 7rh  E d . / . 
1-1 


Similarly,  the  noise  n(t)  can  be  represented  as 
n(t)  = Re{n(t)e^(t)e:ja)ct  > 


where 


n(t)  nc2  (t)  + ng2(t) 


is  the  lowpass  equivalent  envelope  and 


<()  (t)  = tan 


-1  ns(t) 
nc(t) 


is  the  phase.  We  can  then  represent  the  received  signal  as 
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Figure  1.  Optimum  High  SNR  CPMFSK  Receiver 
for  Gaussian  Noise 


r ( t)  = Re 


= Re 


I /2E  j4>  (tl  jto  t 

[t/“^  * di' V ♦ " ««•  }c 

J ja(t)  ju  tl  f ju.O 

|^v(t)  e e j = Re  Sr(t)  e r 


where 

v(t)  is  the  lowpass  equivalent  envelope  of  the  signal 
plus  noise,  a(t)  is  the  lowpass  equivalent  phase,  and  r(t)  is 
the  complex  lowpass  equivalent  received  signal 

/2Eo  j<Mt) 


r (t)  =y/  -y-  s ( t , dx,  Dk)  + n(t)  e 
s 


Whereas  in  the  bandpass  representation,  the  output  of  one 


of  the  correlators  is 


O <5 


r (t)  s (t,  d , D . ) dt, 
A J 


It  can  quickly  be  shown  that  the  same  output  can  be  represented 
using  the  baseband  equivalent  signals  as 


= 1/2  Re  S / r(t)  s*  (t 


, dx.  D.)  dt| 


where  s*  (t,  d, , D.)  indicates  the  complex  conjugate  of 

A D 

s(t,  Dj). 

The  alternate  complex  baseband  representation  of  the 
received  signal  as  envelope  and  phase  as  indicated  in  Eq.  9 is 

ja  (t) 

r (t)  = v (t)  e . (i: 

This  representation  is  useful  when  examining  receiver  effects. 
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2.2  Receiver  Representations 

The  receiver  structures  which  were  investigated  during 
this  study  can  be  depicted  as  shown  in  Figure  2.  The  received 
signal  is  processed  in  a prediction  filter  to  limit  the  input 
bandwidth  of  the  additive  noise.  The  predetection  filters 
considered  in  this  report  are  the  physically  realizable  Butter- 
worth,  Chebyshev,  and  Bessel  bandpass  filters.  In  terms  of 
the  complex  baseband  equivalent  signals,  as  shown  in  [8], 

oo 

y(t)  = / h ( x ) r(t-t)  dx  (14) 

— CO 

where  h(t)  is  the  complex  baseband  equivalent  impulse  response 
of  the  filter.  Since  y(t)  is  also  complex,  it  can  be  represented 
in  terms  of  an  envelope  and  phase 

j 6 ( t) 

y(t)  = p(t)  e (15) 

The  zero  memory  nonlinearities  are  included  in  the 
receiver  structure  to  limit  or  delete  large  impulsive  noise 
spikes.  The  output  of  these  nonlinearities  can  similarly  be 
represented  in  envelope  and  phase  form  as 

jC  (t) 

Z(t)  = ^ (t)  e (16) 

where  £(t)  and  C(t)  are  functions  of  the  input  envelope  and 
phase  p(t)  and  B(t).  In  this  study,  we  consider  the  commonly 
used  nonlinearities  of  the  bandpass  clipper  and  the  bandpass 
hole  puncher  or  blanker  as  shown  in  Figure  3.  For  these  non- 
linearities,  the  output  phase  equals  the  input  phase 

C(t)  = B ( t)  (17) 

whereas  the  output  envelope  equals  the  input  envelope  up  to  a 
threshold  and  it  equals  a constant  if  the  input  envelope  exceeds 
that  threshold.  That  is. 
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INPUT  ENVELOPE 
(a)  CLIPPER 


INPUT  ENVELOPE 
(b)  HOLE  PUNCHER 


Figure  3:  Envelope  Transfer  Functions  of  Zero- 
Memory  Bandpass  Nonlinearities 
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p(t) 


p (t) 


c(t) 


1 fch 


p (t) 


> t. 


(18) 


where  is  the  threshold;  K=t^  for  a bandpass  clipper  and  K=0 
for  a bandpass  hole  puncher. 

After  the  nonlinearity,  the  signal  Z(t)  is  correlated 
in  a linear  receiver  with  all  possible  noise-free  signals  pro- 
ducing the  Mn  output  voltages  x. A=l,  . .,  M,  j = l,  ...»  m. 

A 3 

Similarly  to  Eq.  12,  the  output  voltage  can  be  written  as 


= 1/2  Re 


s*  (t,  dA, 


(19) 


The  rest  of  the  processing  on  the  voltages  is  identical  to 

that  which  exists  in  the  linear  receivers  derived  on  the  Gaussian 
noise  channel  [2-4].  In  Section  3 below,  the  results  of  the 
effects  of  the  predetection  filter  on  performance  and  one 
receiver  implementation  based  upon  the  MLSE  algorithm  are 
presented . 
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3.  Receiver  Implementation  Considerations 

There  were  two  separate  phases  of  the  study  concerning 
receiver  implementations.  The  first  was  to  determine  the  effect 
of  preceding  the  receiver  with  a realizable  bandpass  filter.  The 
second  was  to  consider  more  realistic  receiver  implementations 
than  the  mathematical  structure  of  Figure  1.  The  MLSE  algorithm 
receiver  is  discussed. 

3.1  Predetection  Filter  Effects 

Three  generic  types  of  realizable  filters  were  considered: 
the  Butterworth,  Chebyshev,  and  Bessel.  In  addition,  filter 
order  of  2,  3 and  4 were  evaluated  for  each  filter  type.  For  the 
Chebyshev  filter,  a .1  dB  passband  ripple  was  considered. 

The  transfer  functions  and  impulse  responses  of  these 
filters  can  be  found  in  a number  of  references,  however  they  will 
be  included  here  for  completeness.  Table  1 gives  the  transfer 
functions  and  impulse  responses  of  Butterworth  filters  normalized 
to  a lowass  3 dB  bandwidth  of  1 rad/sec.  Similarly,  Table  2 
gives  the  same  information  for  Chebyshev  filters  with  passband 
ripple  of  .1  dB.  For  Bessel  filters,  the  normalized  lowpass 
transfer  functions  are  not  normally  designed  for  a 3 dB  bandwidth 
of  1 rad/sec.,  but  rather  for  a constant  group  path  delay.  Cor- 
respondingly, the  transfer  functions  and  impulse  responses  must 
be  frequency  scaled.  It  is  necessary  to  know  the  normalized  3. 
dB  bandwidths.  In  Table  3 is  shown  the  normalized  transfer  func- 
tion, the  frequency  transformation  b,  and  the  corrected  impulse 
responses.  The  Laplace  transform  relationship 

H (bs)  ► i h (^)  (20) 


is  used. 
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Table  1:  Transfer  Functions  and  Impulse  Responses 
for  Lowpass  Butterworth  Filters 


i 


Order 

Transfer  Function  H(s) 

Impulse  Response  h(t) 

2 

1 

s2  + 1.414s  + 1 

1 . 414e~ ' 707tsin  ,707t 

3 

1 

e t+1.155e  2tsin 

(s  + 1) (s2  + 1.848s 

+ 1) 

( . 866t+l .047) 

1 

0 - . 924t  • 

2.416e  sin 

A 

(. 

•3Q  9+-_i_  O Q 9 \ 

(s2  + 

JO  L > • J j 

. 765s+l) (s2  + 1.848S+1) 

-e 

382ts in (. 924 t+ 1.178) 

Table  2:  Transfer  Functions  and  Impulse  Responses 

for  Lowpass  Chebyshev  Filters  (.1  dB  ripple) 


Order  Transfer  Function 

H(s)  Impulse  Response  h(t) 

0 3.27 

s2  + 2.36s  + 3.31 

2.37e-1*18tsin(1.38t) 

3 1.59 

.941e~’989t+1.014e~"485t 
sin ( 1 . 206t-l . 19 ) 

(s  + .969)  ( s2+ . 97s  + 1.69) 

4 -817 

, . - . 638t 

1 . 426e 

sin ( . 464t  + .285) 

(s2+. 528s+l. 33) (s2+l. 275s+.622 

) - . 264t 

- . 5 9e 

sin ( 1 . 123t  + .747) 

In  choosing  a predetection  filter,  two  requirements  are 
desired.  First  it  is  desirable  in  many  applications  that  the 
filter  bandwidth  be  narrow  so  that  the  noise  power  which  reached 
the  correlators  is  small.  If  zero-memory  nonlinearities  are  used, 
it  has  been  found  that  wider  bandwidths  are  more  desirable. 
Secondly,  the  filter  bandwidth  should  be  large  enough  to  effec- 
tively pass  the  signal  without  distortion.  The  block  diagram  of 
Figure  4 depicts  the  technique  which  was  used  in  this  study  to 
determine  one  effect  of  the  predetection  filter. 

One  of  the  Mn  complex  baseband  equivalent  signals 
s(t,d^,  D^)  is  chosen  randomly.  It  is  passed  through  a filter 
with  impulse  response  h(t).  The  output  signal  is  then  correlated 
with  an  unfiltered  version  of  itself.  The  delay  compensation  is 
necessary  to  correct  for  the  group  path  delay  through  the  filter. 
From  [9],  the  group  path  delays  for  the  nine  filters  of  considera- 
tion is  shown  in  Table  4.  These  group  path  delays  are  for  the 
normalized  bandwidths  of  1 rad/sec.  For  other  bandwidths,  they 
must  be  adjusted  by 


t 


new 


told 


old 

B 

new 


(21) 


where  t ^ is  the  group  path  delay  as  indicated  in  Table  4,  and 
Bnew  is  the  new  bandwidth  in  rad/sec.  B0^  = 1 rad/sec. 

Clearly,  if  there  were  no  filter  in  Figure  4,  the  output 
of  the  correlation  would  be 


nT 

wlk  “ 2 Re<j^  f S|s(t,dlf  Dr)|2  dt 
Hence,  we  define  a correlation  degradation  as 


nT, 


(22) 
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Figure  4.  Technique  to  Determine  Effect  of  Predection  Filter 


Table  4.  Group  Path  Delays  (sec)  for  Filters 

Normalized  to  Bandwidths  of  1 rad/sec. 


Filter  type 

Order 

2 

3 

4 

Butterworth 

1.6 

2.3 

3.1 

Chebyshev 

.8 

1.85 

2.8 

Bessel 

1.45 

1.85 

2.19 

Ik 

Cd  = 20  log  (23) 


It  is  difficult  to  directly  relate  to  the  error  rate 
performance  because  noise  is  not  included  in  Figure  4 and  the 
effect  on  other  correlators  (refer  to  Figure  1)  is  not  being  in- 
cluded. Nevertheless,  can  be  used  as  a guide  in  choosing  pre- 
detection filters  and  bandwidths  such  that  negligible  distortion 
occurs  in  the  signal.  Figures  5,6,  and  7 show  as  a function  of 
3 dB  bandwidth.  From  the  figures,  it  can  be  seen  that  is 

degraded  by  less  than  1 dB  if  the  one-sided  3 dB  bandwidth  is 
greater  than  .6/Tg.  Furthermore,  the  effect  of  different  filters 
appears  to  be  a minimal  one.  That  is,  as  long  as  the  filter  band- 
width is  moderately  wide  (> ,6/T^) , it  does  not  seem  to  matter  too 
much  what  filter  type  or  filter  order  is  being  used. 

For  quaternary  CPMFSK  of  large  deviation  ratios  (i.e., 
h = .6  or  h = .8),  the  transmitted  spectra  is  considerably  larger 
than  for  h = .4  (see  Section  6.)  Thus,  the  3 dB  bandwidths  of 
predetection  filters  should  be  correspondingly  larger  in  these 
cases  to  ensure  that  minimal  distortion  and  degradation  is  intro- 
duced by  the  filter. 
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BINARY  CPFSK 
h =2/3 


EQUIVALENT  LOWPASS  BANDWIDTH  (3db) 

4357  77E 

Figure  5.  vs.  3 dB  Bandwidth  for  Binary  CPMFSK  and 

Various  Filters:  (a)  Butterworth;  (b)  Chebyshev 
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Figure  6. 
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C,  vs.  3 dB  Bandwidth  for: 
d 

(a)  Binary  CPFSK  with  Bessel  Filter 

(b)  Quaternary  CPMFSK  with  Butterworth  Filter 
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QUATERNARY  CPMFSK 
h = .4 

(a)  CHEBYSHEV  FILTERS 

O 2nd  ORDER 
A 3r^  ORDER 
□ 4th  ORDER 
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QUATERNARY  CPMFSK 
h = .4 


(b)  BESSEL  FILTERS 
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EQUIVALENT  LOWPASS  BANDWIDTH  (3dB) 
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Figure  7.  vs . 3 dB  Bandwidth  for  Quaternary  CPMFSK  and 

Various  Filters  (a)  Chebyshev;  (b)  Bessel 


3.2  Use  of  the  MLSE  Algorithm 


It  was  required  in  the  study  to  briefly  investigate  more 
realistic  receiver  implementations.  One  attractive  alternative 
is  the  MLSE  algorithm,  or  as  it  is  sometimes  more  commonly  known, 
the  Viterbi  Algorithm 


As  indicated  in  Eq.  1, 


cos 


0)  t + 
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dj^  iTh  (t-  (i-1)  Ts 
s 


i-1 

2 

j=o 


Trhd  . 
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(24) 


where  the  definition  of  terms  is  given  in  Section  2.  d^  is  the 
digital  modulation  which  is  +1 , +3 , . . . ,+ (M-l)  for  M-ary  signaling. 
Figure  8 shows  the  phase  trajectory  for  quaternary  signaling 
(d^  = +l,+3).  If  h is  a rational  number,  the  tree-like  depiction 
of  Figure  8 will  fold  upon  itself  modulo  2tt  to  give  a trellis- 
like depiction.  Figure  9a  shows  the  trellis  depiction  for  binary 
signaling  and  h = 2/3.  Figure  9b  shows  the  trellis  depiction  for 
quaternary  signaling  and  h = 2/5.  Reference  [11]  indicates  these 
respective  values  of  h give  good  performance  for  Gaussian  noise 
and  tight  spectral  occupancy  characteristics.  The  values  in 
parentheses  indicate  the  digital  d^  = +1,  +3  which  causes  that 
particular  path  in  the  trellis.  It  may  be  easier  to  depict  the 
phase  trajectories  as  a phase  state  space.  Figures  10a  and  10b 
respectively  show  that  point-of-view  for  the  binary  and  quater- 
nary modulation. 

We  wish  to  find  the  path  through  the  trellis  or  state 
space  which  most  closely  represents  the  received  signal.  The 
details  of  the  MLSE  or  Viterbi  Algorithm  will  not  be  given  here. 
The  interested  reader  is  directed  toward  [7,  10]. 

For  N nodes  in  the  state  space  diagram,  there  will  be  N 
path  metrics  and  N path  histories.  In  addition  there  are  M*N 
branch  metrics  which  are  basically  the  correlations  of  the 
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b.  Quaternary  Trellis  with  h = 2/5 
Figure  9:  Phase  Trellises  for  Some  CPMFSK  Signals 
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a.  Binary  with  h = 2/3 


Figure  10:  Phase  State  Representation  of  Phase  Trajectories 
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received  signal  with  the  M branches  emanating  from  each  of  the 
N nodes.  If  we  denote  these  branch  metrics  as 

T 

bmn  = f S r(t)  cos  [wct  + t + n]  dt  (25) 

• s 

o 

where  m = - (M-l),  - (M  - 3) , . . . ,1, 3, . . . ,M  - 1 (M  terms)  and  n 
represents  one  of  the  N phase  nodes. 

The  updating  of  the  path  histories  can  be  expressed  as 

new  P = max  [old  P + bmn^]  (26) 

n'  n" 

where  new  Pn  is  the  new  updated  path  metric  for  the  n'"  phase 

node,  the  maximization  is  accomplished  over  only  those  n'  which 

have  paths  from  n'  to  n and  b ^ is  the  branch  metric  from  n'  to 

mn 

n.  Old  Pn>  of  course  denotes  the  old  path  metric  at  n'.  Referring 
to  Figures  9 and  10,  we  can  see  that  n'  refers  to  all  of  the  other 
n - 1 phase  nodes  because  there  are  paths  from  all  other  phase 
nodes  to  n for  any  n. 

To  perhaps  more  clearly  explain  Equation  26,  for  the 
binary  example  of  Figures  9a  and  10a,  the  three  updates  of  path 
metrics  are 

new  Pq  = max  [(old  P2/3tt  + b-i,2/3v)  ' (old  P-2/3tt+  bl,-2/3TT)1 

neW  P2/3tt  = max  [(old  P0  + bl,0>'  (old  P-2/3tt  + b-l,-2/3Tr)  ] (27) 

new  P-2/3tt  = max  [(old  P2/3tt+  bl,2/37T)  ' (old  P0  + b-l,0)] 

For  each  P , after  the  n'  which  accomplishes  Equation  26 
is  found,  the  path  histories  are  updated.  The  path  history  update 
can  be  expressed  as 
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d 


(28) 


new  H = old  H ' , m 
n n 

which  shows  that  the  new  path  History  H is  produced  by  appending 
to  the  old  Hn^  the  data  digit  m which  corresponds  to  the  branch 
from  n'  to  n.  The  n'  used  is  of  course  the  one  which  satisfies 
Equation  26  for  each  n. 

Ideally,  we  would  like  to  collect  path  metrics  and  path 
histories  for  the  entire  duration  of  a message.  In  practice, 
it  has  been  found  that  it  is  not  necessary  to  wait  until  the  end 
of  a message  to  start  making  symbol  estimates,  but  after  some 
delay  D estimates  can  be  made.  Hence,  once  the  length  of  the 
histories  equals  D,  a path  is  chosen  which  has  the  largest  path 
metric  and  the  oldest  symbol  in  that  path  history  is  chosen  as 
the  symbol  estimate  of  a previous  received  symbol  with  delay  D. 

The  choice  of  D depends  upon  the  application,  and  was 
investigated  by  simulation. 

Monte  Carlo  simulations  programs  were  developed  and 
exercised  for  the  binary  case  with  h = 2/3,  corresponding  to 
Figures  9a  and  10a.  Gaussian  noise  was  used,  and  Figure  11 
shows  the  simulated  probability  of  bit  error  results,  as  circles 
on  the  graph.  These  results  are  for  a Viterbi  Algorithm  delay 
D = 3.  Also  shown  is  the  theoretical  bound  from  (11]  which 
approximates  the  error  performance  of  a correlation  receiver 
which  observes  three  bits  (n  = 3)  before  a decision  is  made. 

As  can  be  seen,  the  performance  of  both  techniques  are  equivalent. 

If  the  Viterbi  Algorithm  delay  D is  increased,  perfor- 
mance will  improve,  and  if  decisions  are  delayed  until  the  entire 
message  has  been  sent,  the  best  performance  possible  will  result. 
However,  this  is  not  truly  necessary  for  eventually  a point  is 
reached  for  which  an  increasing  of  D beyond  that  point  achieves 
virtually  the  same  results. 
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O SIMULATION  POINTS 
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Figure  11:  Probability  of  Bit  Error  for  Viterbi  Algorithm 
Demodulation  - Binary  Signaling  h = 2/3 


Figure  12  shows  some  simulated  results  for  D = 2,  3,  and 
5.  Results  were  also  obtained  for  D = 7,  but  they  are  graphically 
indistinguishable  from  the  D = 5 points.  This  same  improvement 
can  be  achieved  by  using  a correlation  receiver  which  observes 
five  bits  (n  = 5) , but  this  requires  an  increase  in  equipment 
complexity,  whereas  the  Viterbi  Algorithm  Receiver  does  not. 

The  Viterbi  Algorithm  Receiver  is  only  practical  when 
the  deviation  ratio  h is  a rational  number  such  that  the  state 
space  depiction  of  Figure  10  has  a small  number  of  phase  nodes. 
Fortunately,  CPMFSK  is  a modulation  which  is  robust  enough  such 
that  these  "convenient"  values  of  h result  in  good  performance. 

Equation  25  implies  that  M*N  correlators  are  required. 
This,  in  fact,  is  unnecessary  as  we  can  see  by  rewriting  Equation 
25  as 


b = cos  (n)  / 
ran  n 


r(t)  cos  [u)  t + 
c 


-t]  dt 


- sin  (n)  S 


r ( t)  sin  [cj  t + 
c 


-t]  dt 


Since  the  N phase  nodes  are  known,  it  can  be  seen  that 
only  two  quadrature  correlators  are  required  for  each  m.  Denoting 
these  by 


x = f 
m 

o 


r(t)  cos  [ u>  t + 
c 


-t]  dt 


r s • r . , mTTh.  , 

ym  = / r ( t)  sin  [w  t + t] 

o s 


b„  can  be  written  as 
mn 


b = cos  (n)  * x - sin  (n)  • ym 

mn  m 


J 


Hence,  a total  of  2M  correlators  are  required  for  M-ary  modulation 
rather  than  N-M.  Figure  13  shows  the  functional  block  diagram  of 
a binary  receiver  ( m = +1)  with  h = 2/3  (n  = 0 , | tt  , - | it)  . 
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Figure  13:  Functional  Block  Diagram  of  Binary  Receiver  with 
Viterbi  Algorithm  Demodulation 


4.  Noise  Models 


Three  noise  models  were  used  during  the  course  of  the 
study:  two  are  suitable  for  theoretical  analysis,  and  one  for 

simulation  purposes.  Two  techniques  were  used  to  specify  the 
models.  First,  the  models  were  required  to  have  finite  Vd  ratios, 
which  is  defined  in  terms  of  the  rms  to  average  envelope  levels 
as 

/ir2 

vd  = 20  log  <32> 

Vd  is  a commonly  used  parameter  to  specify  the  degree  of  impul- 
siveness of  the  noise.  In  Eq.  32,  p_.  is  the  jth  moment  of  the 
envelope  of  the  noise.  In  terms  of  Eqs.  6 and  7, 

oo  -i 

Mj  = fQ  nJ  pn  (n)  dn  (33) 

where  PN(n)  is  the  envelope  probability  density  function  (pdf). 

For  Gaussian  (Rayleigh  envelope)  noise,  Vd  is  1.05  dB,  whereas 
common  values  are  2-4  dB  for  HF  channels  and  10-15  dB  for  VLF 
channels . 

The  second  technique  used  to  specify  the  models  is  the 
empirical  APD  curves  of  [6].  Models  which  match  these  curves 
are  preferred  since  they  are  obtained  from  measured  data.  The 
APD  curve  is  a plot  of  the  probability  that  the  noise  envelope 
exceeds  a threshold  R,  i.e. 

Pr(n>R)  = C PN(n)  dn  (34) 


* 
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4.1  Analytical  Noise  Models 

In  this  study,  two  basic  noise  models  which  achieve  the 
joint  goals  of  analytical  tractability  and  correspondence  with 
empirical  APD  curves  and  ratios  were  considered.  One  of  these 
models  has  been  previously  proposed  by  Hall  [12] , although  for 
some  Vd  ratios,  modification  of  the  model  was  required.  The 
other  model,  called  a mixture  model,  represents  the  envelope 
noise  process  as  a mixture  of  a Rayleigh  process  and  an  impulsive 
process  to  be  defined. 

4.1.1  The  Hall  Model 


The  Hall  model  can  be  represented  by  its  envelope  pdf 

(35) 


P(n)  - 


[n2+v2] 


0 + 1 


where  9 is  an  integer  ]>  2 which  is  selected  to  match  first  order 
statistics.  y is  a parameter  related  to  the  variance  of  the 

. l- 

noise  if  this  variance  exists.  The  j moment  of  this  process 
does  not  exist  if  j > 0-1  so  that  in  these  cases  modifications  of 
the  process  is  required  to  define  a Vd  ratio.  Hall  suggests  one 
method  of  modifying  the  model  so  that  finite  moments  result. 

The  method  proposed  here  is  described  below  and  differs  from 
Hall's  suggestion.  Marginal  density  functions  which  represent 
the  inphase  and  quadrature  voltages  nc  and  ng  can  be  derived. 

For  this  model  APD ' s can  easily  be  determined  and  as  seen  in 
Figures  14  and  15  the  match  between  the  theoretical  and  empirical 
curves  is  quite  evident. 

The  two  parameters,  0 and  Y,  of  the  Hall  model  whose 
envelope  pdf  is  given  in  Eq.  35  above  can  be  selected  to  provide 
good  fits  to  the  measured  distributions  of  atmospheric  noise. 
Integral  values  of  0 in  the  range  2 <_  0 £ 5 are  commonly  used. 

The  distributions  for  which  0=4  and  5 represent  moderately 
impulsive  noise  with  = 3 and  2 dB  respectively  and  are  parti- 
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cularly  suited  for  modeling  the  first  order  noise  statistics  in 
mid-ocean  regions.  However,  the  moments,  of  the  Hall  distri- 

butions diverge  for  j:>0-l.  Thus,  the  second  moment  of  the  dis- 
tribution for  0=3  and  the  first  and  second  moments  for  the 
0=2  distribution  do  not  exist.  Consequently,  for  these  dis- 
tributions, the  ratio  as  defined  in  Eq.  32  is  undefined  and 

the  distributions  cannot  be  constrained  to  match  the  V,  ratio  of 

d 

the  empirical  data.  However,  truncation  and  appropriate  normali- 
zation of  these  distributions  force  the  moments  to  converge  and 

result  in  a set  of  distributions  whose  V,  ratio  is  defined  in 

d 

terms  of  0 and  the  truncation  level. 


The  truncation  and  normalization  is  accomplished  by 
defining  an  envelope  pdf  as 


PE(n)  = 


b(0-l)Ya  1n 


- o < n < e 


r 2 ^ 2 

|_n  + y 


(0+l)/2 


and  selecting  b to  provide 


/EPp(n)  dr)  = 1 
0 E 


(36) 


(37) 


These  truncated  distributions  with  0=2  and  3 and  the  unmodified 
distributions  with  0=4  and  5 constitute  a family  of  curves 
which  can  be  used  to  represent  a wide  range  of  measured  data  by 
exhibiting  the  same  ratio  as  the  data  and  providing  good  fits 
to  the  observed  APD’s.  The  family  is  well-defined  analytically 
and  expressions  can  be  obtained  in  closed  form  for  the  moments  of 
the  envelope,  APD's  and  marginal  pdf's  of  the  quadrature  voltage 
components.  These  statistics  are  summarized  in  Table  5. 
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Table  5.  First  Order  Statistics  of  Truncated  Hall  Model 
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The  marginal  pdf  is  obtained  from  the  envelope  pdf 
through  the  transformation  nc  = ti  cos  <j>  and  ng  = n sin  <p 
where  4>  is  independent  of  n and  is  assumed  to  be  uniformly 
distributed  in  the  interval  (0,  2it).  The  details  of  the  deriva- 
tion of  the  expressions  in  Table  5 are  provided  in  Appendix  B. 

The  correspondence  between  the  theoretical  APD  express- 
ions which  are  presented  in  Appendix  B and  the  empirical  CCIR 
results  given  in  [6]  are  graphically  displayed  in  Figures  14 
and  15.  A normalized  parameter  A is  defined  as 

A = 20  log  (-2 — ) (38) 

nrms 

where  nrms  = and  the  figures  display  the  probability  that  A 

is  exceeded.  These  parameters  correspond  to  those  used  in 
reference  [6],  and  the  empirical  results  displayed  in  [6]  are 
shown  as  isolated  points  in  the  figures.  The  solid  or  dotted 
lines  correspond  to  the  theoretical  APD  expressions  as  stated 
in  Table  5. 

Figures  14  and  15  show  the  results  for  the  Hall  and 

truncated  Hall  model.  Only  for  = 2 dB,  for  which  6=5,  can 

the  original  Hall  model  be  used.  For  all  other  results  in  the 

figures,  it  was  necessary  to  use  the  truncated  models  with 

either  0 = 3 or  6 = 2.  It  was  found  experimentally  that  8=3 

gives  the  best  fit  to  the  CCIR  curves  when  4 < V,  < 6 . For  V, 

— d — d 

larger  than  6,  0 = 2 gives  the  best  fit.  Clearly  the  fit  is 
good  for  all  the  curves. 

This  noise  model  has  also  been  used  in  a study  to 
evaluate  low  probability  of  intercept  communications  in  a broad- 
band spread  spectrum  system  using  PSK  as  a modulation  [13]. 
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100  X PROBABILITY  A IS  EXCEEDED 


Figure  15:  APD  for  Empirical  and  Truncated  Hall 
Model  (V(j=10,  12,  and  14) 
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4.1.2  The  Mixture  Model 


The  mixture  model  considered  here  is  represented  by  its 
envelope  process 

n (t)  = -■(--)-2-+-z  (t)  + s(t)  ~ z(t)  u ( t ) (39) 

where  s(t)  is  an  impulsive  process,  z(t)  is  a Rayleigh  process, 
and  u ( t ) is  a switching  process  which  is  independent  of  s(t)  and 
z(t)  and  is  +1  with  probability  e and  -1  with  probability  1-e. 

In  other  words,  when  u(t)  = 1,  the  envelope  is  represented  by  an 
impulsive  process  and  when  u(t)  = -1  it  is  represented  by  a 
Rayleigh.  This  type  of  process  has  been  suggested  before  [14,15], 
although  the  impulsive  part  of  the  process,  s(t),  in  this  study 
allows  a greater  amount  of  closed  form  analysis. 


The  processes  z(t)  and  s(t)  are  assumed  to  be  quasi- 
stationary and  can  be  represented  by  their  pdf's 

2 


P(z)  = 


exp 


2a 


z > 0 


(40) 


p (s)  = 


(ka) 


r+1 


F (r)  2 


r-1 


K 


1-r  [ka]' 


r > 0 


(41) 


2 

where  z is  a Rayleigh  variable,  2a  is  the  second  moment  of  z, 
k controls  the  ratio  of  power  between  the  impulsive  and  Rayleigh 
portions  of  the  noise,  and  r is  related  to  the  impulsiveness  of 
s.  T(r)  is  the  Gamma  function  and  K^_r(  ) is  a modified  Bessel 
function  of  order  1-r.  For  integer  r this  form  of  impulsive  func- 
tion has  been  suggested  previously  [16].  If  r = 1,  the  envelope 
density  function  of  s corresponds  to  identically  distributed 
quadrature  variables  with  a Laplace  pdf;  hence  s will  be  called  a 
generalized  Laplace  variable.  It  is  surprisingly  rich  in  analy- 
tical usefulness,  and  most  of  the  results  of  interest  to  communi- 
cations engineers  can  be  expressed  in  closed  form. 


r 


The  moment  of  s can  be  found  from  equation  (41)  using 

item  6.561-16  in  [17],  i.e., 

j 


E{sj}  = (2y^y  rd+i)  r (r+j) 


(42) 


Specifically  the  mean  is 


E { s } = 


and  the  second  moment  is 


ko  (t+^-) 
T (r) 


(43) 


E{s2}  = 4r (ko ) 2 


The  probability  that  s exceeds  some  fixed  threshold  R is 


(44) 


Pr{s  > R]  = f p (s ) ds 
R 


- (J)  r 


K 


yko  J 2r-1r (r)  * rVk°/  <45> 


If  it  is  assumed  that  s is  the  envelope  of  two  quadrature  varia- 
bles x and  y,  the  pdf  of  x (or  y)  is 


p (x)  = 


x 


H 


r+\  . 

(ko)  ^ /it F ( r ) 2 


r Ki 

r — — — 

^ 2 


r(^) 


(46) 


Eqs.  (45)  and  (46)  are  derived  in  Appendix  C.  Returning  to  the 

t hi 

mixture  model,  it  is  clear  that  the  j moment  of  n is 


u.j  = E{n^|u  = 1}  Pr{u  = 1)  + E(n^|u  = -1}  Pr{u  = -1} 


= e E{  s-* } + ( 1-e  ) E { z ^ } 


(47) 


Specifically,  the  mean  is 
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Similarly,  the  APD  is 

Pr(n  >R)  = Pr{n  > R|u  = 1}  Pr{u  = 1}  + Pr{ri  > r|u  = -1} 

Pr{u  = -1}  (51) 


Figures  16  and  17  display  the  results  for  the  mixture 
model.  For  = 2 , this  model  does  not  achieve  as  good  a fit 
as  does  the  Hall  model.  Nevertheless  there  is  still  a clear 
correspondence  between  the  theoretical  APD  curves  as  derived  in 
the  model  and  the  empirically  determined  results  from  [6]. 

There  are  three  parameters,  namely  r,  e,  and  k,  to  pick 
in  the  mixture  model.  The  parameters  were  chosen  to  best  match 
the  empirical  results  with  the  constraint  that  the  theoretical 
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Vd  ratio  as  defined  in  Eq.  (50)  be  met.  It  can  be  seen  that 
as  the  noise  becomes  more  impulsive,  i.e.,  the  Vd  ratio  get 
larger,  r becomes  smaller  and  e and  k get  larger. 

It  is  of  interest  to  note  that  all  the  moments  of  x 
and  the  exceedance  probabilities  of  x can  also  be  derived  where 
x is  a quadrature  variable  of  the  impulsive  envelope  s(Eq.  46). 

Namely,  the  moments  are  seen  to  be 


E{xj}  = /"  xj  p ( x ) dx  =< 


0 , j odd 

( 2k o ) 

77T 


TrT  * r(^T^ ) r(r+2)  ' 


(52; 

j even 


Regarding  the  exceedance  probabilities, 


/ p ( x ) dx  = / 
R R 


xr'l 


(ko ) r+2  /it  T (r)  2r  2 7 


T Ki 

r-~  ^-r 


(— ) 
vkcr 


R 


r 2 


/ 

0 


r+i  r -i 

(ko)  z /n  r ( r ) 2 ' 


(fc) 


(53) 


From  [17,  item  6.561(4)]  with  a suitable  change  of  variable,  this 
results  in 


(54) 


where  ( * ) is  the  Modified  Struve  Function. 
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4.2  Simulation  Noise  Model 


The  simulation  noise  model  used  in  this  study  is  based 
upon  the  Crichlow  graphical  model  for  the  APD  of  the  noise 
envelope  [18 ] . 

The  Crichlow  model  works  with  a special  type  of  probabil- 
ity graph  paper;  paper  on  which  the  power  Rayleigh  functions 
plot  as  straight  lines  (a  Rayleigh  function  plots  as  a straight 
line  with  slope  = -0.5).  The  coordinate  transformations  for  the 
probability  paper  are  (see  Appendix  D) : 


x'  = -20  log  (-In  x) 
y'  = 20  log1Q  y 


(55) 


where  x = the  precentage  of  time  for  which  the  ordinate  is 
exceeded;  and  y=n/n 

rms 

On  this  probability  paper,  the  APD  of  atmospheric  noise 
can  be  represented  by  a three-section  curve  (see  Figure  18) . The 
lower  region  of  the  curve,  representing  random  low  voltages  and 
high  probabilities,  approaches  a Rayleigh  distribution.  It  can 
be,  therefore,  approximated  by  a straight  line  (R) . The  higher 
region  of  the  curve,  representative  of  impulsive  high  voltages 
and  low  probabilities,  approaches  a power  Rayleigh  distribution. 

It  can  also  be  approximated  by  a straight  line  (PR) . The  center 
region  of  the  curve  corresponds  to  an  arc  of  a circle  tangent 
to  the  two  straight  lines. 

The  arc  is  also  tangent  to  the  line  (T) , which  is  parallel 
to  line  (BI)  which  bisects  the  acute  angle  formed  by  the  inter- 
section of  the  Rayleigh  and  power  Rayleigh  lines. 

Four  parameters  are  necessary  in  order  to  specify  a unique 
pair  of  lines  and  an  arc.  They  are: 
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Figure  18.  Crichlow  Graphical  Model  of  APD 
of  Atmospheric  Noise 


1)  slope  of  the  power  Rayleigh  line; 

2)  point  through  which  the  power  Rayleigh  line  passes; 

3)  point  through  which  the  Rayleigh  line  passes  (the 
Slope  is  known  to  be  -1/2) ; 

4)  parameter  determining  the  radius  of  the  circular 
arc . 

Crichlow  defined  the  four  parameters  as  follows 

1)  X = -2s,  where  s is  the  slope  of  the  power  Rayleigh 
line ; 

2)  C(dB)  = the  dB  difference  between  the  power 
Rayleigh  line  and  the  Rayleigh  line  at  p = 0.01; 

3)  A(dB)  = the  dB  value  of  the  Rayleigh  line  at 
p = 0.5; 

4)  B(dB)  = the  dB  difference  between  the  y'-axis 
intercepts  of  lines  (BI)  and  (T) . 

Experimentally  measured  APD's  indicate  that  parameter  B 
is  linearly  related  to  first  order  to  the  parameter  X by 

B = 1.5  (X-l)  (56) 

Thus,  once  the  values  of  parameters  X,  C and  A are  known, 
a unique  APD  can  be  constructed. 

The  X,  C and  A parameters  vary  according  to  the  value 
of  the  APD.  Wilson  [19]  has  calculated  the  values  of  these 
parameters  for  values  of  4.0  through  30.0.  Parameter  values 
for  = 2.0  and  3.0  were  determined  by  extrapolating  upon 
Wilson's  values  along  with  experimental  curve  fitting. 

Two  computer  subroutines  accomplish  the  simulation  of 
the  random  atmospheric  radio  noise.  They  are  called  "CONSTR"  and 
"ANOISE" . CONSTR  constructs  the  APD  for  a desired  value, 
while  ANOISE  uses  the  parameters  computed  in  CONSTR  to  generate  a 
random  number  with  APD  of  atmospheric  noise.  The  simulation  is 
broken  up  into  two  separate  procedures  since  in  order  to  simulate 
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atmospheric  radio  noise  with  a certain  Vd  value  subroutine  ANOISE 
will  be  called  repetitively;  while  the  parameters  for  the  APD 
for  this  certain  value  have  to  be  calculated  only  once. 

Subroutine  CONSTR  accepts  for  input  the  value  of 
variable.  It  then  looks-up  the  X,  C and  A parameters  corres- 
ponding to  the  value.  At  present,  the  subroutine  has  memory 
of  the  three  parameters  for  values  of  2.0  through  16.0,  in 
steps  of  1.0.  Parameters  for  the  intermediate  values  of 
are  computed  by  linear  interpolation.  Using  the  three  parameters, 
CONSTR  then  "constructs"  the  APD  according  to  the  Crichlow  graph- 
ical model. 

The  output  of  the  CONSTR  subroutine  consists  of  the  two 
parameters  defining  the  borders  between  the  three  regions  of 
the  APD;  and  the  seven  parameters  necessary  for  the  construction 
of  different  sections  of  the  APD: 

1-2)  slope  and  y' -intercept  of  the  Rayleigh  line; 

3-4)  slope  and  y' -intercept  of  the  power  Rayleigh  line; 

5-7)  x'  and  y'  coordinates  of  the  center  of  the  circular 
arc,  and  its  radius. 

The  nine  parameters  computed  by  the  subroutine  CONSTR 
constitute  the  input  to  the  subroutine  ANOISE. 

In  subroutine  ANOISE,  first  a uniform  random  number 
between  0 and  1.0  is  generated;  this  number  represents  the 
abscissa  value.  It  is  then  transformed  into  the  x'  coordinate. 
Decision  is  made  as  to  within  which  of  the  three  regions  of  the 
APD  the  random  value  falls,  and  the  y'  coordinate  (n/nrms)  is 
computed  using  the  appropriate  equation.  The  value  of  n/hrms 
is  then  transformed  from  dB  to  volts.  If  an  instantaneous  quadra- 
ture voltage  is  desired,  the  resulting  n/nrnis  voltage  is  multiplied 
by  cos  2~u  where  u is  another  random  number  between  0 and  1. 
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The  above  pair  of  subroutines  was  thoroughly  tested  with 
a total  of  20,000  random  numbers  generated  by  CONSTR  and  ANOISE, 
for  each  of  ten  different  Vd  values.  The  APD's  were  then 
determined  for  the  ten  types  of  noise.  These  APD's  of  atmos- 
pheric radio  noise  as  determined  by  experimental  measurements  [6]. 
This  can  be  seen  in  Figure  19  for  V ratios  of  2,  6,  10,  and  14 
dB. 
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5.  Error-Rate  Results 

This  section  delineates  the  error-rate  results  of  the 
study.  It  was  desired  to  have  both  theoretical  and  simulated 
estimates  of  the  error  rate,  but  as  mentioned  in  the  introduc- 
tion, budgetary  and  time  constraints  precluded  the  complete 
development  of  the  theory.  However,  the  simulation  phase  of 
the  study  was  successful  and  does  indicate  that  CPMFSK  can 
achieve  lower  error  rates  than  MSK  when  both  are  being  processed 
at  the  same  signal- to-noise  ratio.  Subsection  5.1  below  out- 
lines the  theoretical  considerations  up  to  their  present  stage 
of  development.  Much  of  the  detail  calculations  are  consigned 
to  Appendix  A.  Subsection  5.2  presents  the  bulk  of  the  simula- 
tion work  and  results. 


5.1  Theoretical  Error  Rate  Considerations 

The  theoretical  error  rate  determinations  start  with 
the  union  bound  approach  [2,4].  Referring  to  Figure  1,  and 
given  d^  = v and  {d2,...,dn)  = D_.  , the  probability  of  an  error 
is  by  the  union  bound 


M 


PrM(e|v,D.)  < Z 


m 

Z 

N=1  J=1 

N^v 


Pr(xNJ  > xvj}‘ 


(57) 


Averaging  over  D_.  results  in 


. M m m 
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and  averaging  over  v gives 


PrM(c)  < R 


1 M 

~ Z 
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m 
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N^v  J=1  j=l  v= 1 


rr(xNJ  >Xvj) 


(59) 
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We  have  now  reduced  the  problem  from  the  determination  of  error 
rates  for  M-ary  communications  to  the  binary  problem  of  the 
probability  that  one  variable  xNJ  exceeds  another  xv j • Referring 
to  Eq.  19,  we  can  see  that  the  binary  probability  can  be 
written  as 

nT  * nT 

Pr  (x  >x  .)  = Pr { Re [ / s z(t)s  (t,d  ,D  ) dt  > / S 
NJ  V]  O N J o 


z (t) s* (t,  dv,  D j ) dt ] } 


(60) 


or 


nT 

Pr(xNJ>xvj)  = Pr { Re [ / S z(t){s*(t,  dN,  Dj) 
-s* (t,  d^ , Dj)  } d t ) } 


(61) 


Eq . 61  indicates  that  the  problem  requires  the  evaluation  of 
stochastic  integrals.  It  is  possible  to  do  this  for  Gaussian 
noise  and  linear  receivers,  however  there  are  few  general 
approaches  that  can  be  used  to  evaluate  these  integrals.  Speci- 
fically, there  are  three  difficulties  with  evaluating  Eq . 61  in 
the  general  case  of  CPMFSK: 


(1) 


If  the  receiver  is  nonlinear,  the 
z(t)  are  difficult  to  determine, 
hand,  if  the  receiver  is  linear, 

ITT 


Z ( t) 


s ( t , dv,  Dj)  + n (t) 


statistics  of 
On  the  other 

e j<Mt) 


and  techniques  can  be  developed  to  estimate  the 
integral.  This  is  outlined  below  in  section  5.1.1. 

(2)  s*(t,  dN,  D_. ) and  s*(t,  d^,  D ^ ) are  functions 

of  time  which  change  appreciably  during  the 
integral  period  nTg.  This  is  different  than  the 
approach  which  could  be  taken  with  PSK  type  mod- 
ulations, for  which  the  baseband  equivalent  repre- 
sentations are  constant  during  symbol  intervals. 
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(3) 


s*(t,  dN,  Dj)  and  s*(t,  dv,  D j ) are  not 
orthogonal  in  general.  In  specific  instances 
where  the  deviation  ratio  h is  a multiple  of  1/2, 
the  different  baseband  equivalent  signals  are  ortho- 
gonal, and  certain  simplifying  calculations  can 
be  made  to  Eq.  61.  This  is  explained  in  Section 
5.1.2. 

5.1.1  Linear  Receiver  Approach 

As  is  well  known  [8]  when  a waveform  z(t)  is  correlated 
as  in  Eq.  61  with  a local  signal  [s*(t,  dN,  D^)  -s*(t,  dv , D_.)J, 
this  is  equivalent  mathematically  to  passing  the  signal  z(t) 
through  a filter  with  impulse  response 


h(t)  = s*(nTs-t,  dN,  Dj)  -s*(nTs  -t,  dy,  D.) 


(62) 


and  sampling  the  output  at  time  nTg.  Furthermore,  if  the 
receiver  is  linear  the  input  can  be  represented  as 


z (t) 


j4>  (t) 

Dj)  + n (t)e 


(63) 


where  n(t)  and  <f>(t)  are  the  envelope  and  phase  of  an  impulsive 
noise  process.  It  is  known  [20]  that  when  an  impulsive  process 
is  passed  through  a filter  whose  noise  bandwidth  is  less  than 
the  spectral  occupancy  of  the  original  process,  the  output  noise 
is  still  impulsive  but  with  a smaller  ratio. 

Hence,  if  we  define  y to  be  the  decision  variable 


nT 


y=R  ef 
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2E 
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S (t,  d. 


Dj)  + n(t)  e 


j<J>  (t) 


s*(t,  dN,  Dj) 


-s* 


(t, 


dt 


(64) 
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where  we  have  suppressed  the  functional  dependence  of  y on  N,  J, 
v,  and  j,  y is  an  impulsive  voltage  function  with  nonzero  mean. 
Thus,  we  can  model  y in  terms  of  the  analytical  noise  models 
of  Section  4.1.  As  can  be  surmised  in  Table  5 and  Eqs.  39  and  46, 
three  parameters  will  suffice  to  model  y;  the  mean  and  variance  of 
the  voltage  and  the  ratio  of  the  envelope.  From  Eq . 64, 


E { y } = 


nT 


^JNJ'  v j-1* 


(65) 


where 


NJ , vj 


1 

nT" 


nT 


Re  fs(t,  d^,  Dj)  s*(t,  d. 


N ' 


V] 


dt 


Similarly,  the  variance  of  y is 


(66) 


Var*y}  = a‘ 


nT 


/ 

o 


S I h ( t)  I2 


(67) 


where  h(t)  is  as  defined  in  Eq . 62,  and  o is  of  course  the 
variance  of  the  input  noise  process.  From  this,  we  can  see  that 


Var { y } = 2on  nTs  (1-  PNJ,vj) 


(68) 


To  determine  the  modeled  ratio  of  y,  we  assume  that  the  input 

ratio  of  n is  known,  and  that  the  RF/IF  input  noise  bandwidth 
is  where 

B1  > fi  . 

S 

The  noise  bandwidth  of  a filter  with  impulse  response  h(t)  is 


r | H (w) I2  dto 


B = 


(69) 


H ( 0 ) 
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where  H(co)  is  the  Fourier  transform  (transfer  function  of  the 
filter)  of  h(t).  Using  Parseval's  Theorem 


_ | h (t) I 

| H (0  ) p 


r ? 

0 | h ( t)  | dt 

nT  _ 

| / S h ( t)  dtp 
o 


(70) 


Again,  it  is  useful  to  note  that  is  a function  of  v,  j,  N,  and 
J,  so  using  the  conversion  formulas  of  [20],  a different 
ratio  will  in  general  be  required  for  each  term  in  Eq . 59. 

Recapping,  for  linear  receivers,  Eq . 59  can  be  rewritten 
as 


. . M m m M 

PFM(e)  < m * m E E E E Pr {y (N, J, v, j ) > 0}  (71) 

FM  Mm  ^ J=1  j=1  V=1 


where  the  mean  of  y(N,J,v,j)  is  given  by  Eq.  65,  the  variance 
by  Eq.  68,  and  the  effective  ratio  by  knowledge  of  Bn  (Eq.  70) 
, and  the  input  ratio. 


This  basic  approach  to  linear  receiver  performance  has 
been  used  before  [5,  21,  22]  for  the  cases  of  PSK,  MSK , and  CFSK 
modulations . 


5.1.2  Nonlinear  Receiver  Approach 

When  zero-memory  nonlinearities  are  used  in  the  receiver, 
the  input  to  the  correlators  (or  matched  filters)  can  no  longer 
be  considered  impulsive.  We  have  chosen  in  this  case  to  model 
the  stochastic  integral  of  Eq . 61  as  a finite  sum  using  the 
trapezoidal  rule. 

nT  r L-l  .. 

y=Re  / s z ( t ) h (nT  -t)  dt  ~ Re j-^( f (0 ) +f  (L)  ) +Tlf  ( <0>( 72 ) 
o s (.z  j,=1  J 
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where  f ( £) 


(73) 


= *<V  h (nTs  -t£) 


(74a) 

(74b) 


and 

and 


1 fc£  _t£-l 


£>1 . 


nT 


L = 


In  Eq.  72,  h(t)  is  given  by  Eq . 62,  and  again  we  have  suppressed 
the  functional  notation  of  y on  v,j,N,  and  J. 


z(t^)  can  be  expressed  in  amplitude  and  phase  notation 
using  Eq.  16  as 


Z(V 


£ ( t£)  e 


j'’  ^fc£* 


(75) 


Similarly,  h(nTg  t^)  can  be  expressed  in  envelope  and  phase  form 
as  (using  Eq.  62) . 


h(nW 


Hence , 


-j6 (tD ) 

s (t^,  d^ , Dj ) - s*(t^,dv,Dj)  = hft^e 

(76) 


, iU(t£}  " <5  (tp ) ] 

f(£)  = C(t„)  h(t„)  e3  * 1 


Using  this  notation  we  can  rewrite  Eq.  72  as 

L L 

y=  t ( t ^ ) cos  (^  (t£) -6  (tj,))  = z a 


1=0 


where  a is  a constant 


£=o 


£ y£ 


a£  < 


2 h ( t^ ) , £=0  and  L 

Th(t£)  , 1 < £ < L-l 


(77) 


(78) 


(79) 
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and 


yz  = C(t^)cos  ( c < ) - 6 ( t£ ) ) 


(80) 


Two  questions  arise  in  the  further  development  of  Eq.  78.  First, 
how  many  samples  or  points  should  be  used  to  approximate  the 
integral,  and  second,  what  are  the  statistics  of  the  samples 
C(t^)  and  C(t^)  and  how  might  they  be  used  to  determine  Pr{y  >0}, 
which  is  the  desired  computation  (see  Eq.  61)  ? All  of  the 
results  in  this  report  assume  that  samples  of  the  signal  at  time 
t^  are  independent  from  samples  at  any  other  time  t , m^j.  This 
implies  that  samples  are  taken  at  the  Nyquist  rate  relative  to 
the  baseband  equivalent  bandwidth  of  the  bandpass  zero  memory 
nonlinearity  and  that  the  noise  is  white.  If  samples  are  taken 
at  any  other  rate,  there  will  in  general  be  correlation  between 
samples  at  t^  and  tm-  This  can  be  handled  in  the  theory,  but  it 
does  complicate  the  development.  (More  is  mentioned  in  this 
regard  in  Appendix  A) . 

If  the  number  of  samples  L becomes  very  large,  the  law  of 
large  numbers  can  be  used  to  argue  that  y in  Eq.  78  is  approxi- 
mately Gaussian.  With  this  assumption,  the  pdf  of  y can  be 
determined  with  knowledge  of  only  the  mean  and  variance  of  y. 
Furthermore,  as  can  be  seen  from  Eq . 78,  knowledge  of  the  first 
and  second  moments  of  each  of  the  samples  y^  would  allow  deter- 
mination of  the  mean  and  variance  of  y.  Namely,  from  Eq.  78 


and 


L 

E (y ) = £ a E{y  } (81) 

SL=o  * *■ 


2 k 2 2 L L 

E{y^}  = E a j E{y.  } + 21  Z a a. 
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(82) 
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Because  of  our  assumption  of  independent  samples,  in  Eq.  81, 

Eiyfy-jt  ~ Ely^jEiy-j}  and 

L 2 

Var(y)  = £ a^  Var(yy)  (83) 

li=o 

If  the  number  of  samples  L is  not  very  large  (L<20) , then 
one  approach  which  can  be  taken  is  to  use  a Gram-Charlier 
expansion  of  the  pdf  of  y.  This  approach  has  been  used  previously 
to  determine  error  rate  estimates  [22,  23]  when  nonlinear  receives 
are  used  on  the  VLF/LF  channel.  This  technique  requires  the 
knowledge  of  higher  moments  of  the  time  samples  y ^ , but  it  is 
a general  approach  which  has  proven  to  be  useful  in  many  diverse 
fields  and  computations  for  which  the  Gaussian  assumption  is 
not  warranted.  Details  of  the  present  development  of  this 
approach  is  given  in  Appendix  A. 

In  general,  a numerical  approach  is  necessary  to  determine 
the  moments  of  the  time  samples  y^.  This  is  further  explained 
in  Appendix  A,  and  it  seems  to  be  a result  of  the  following: 

The  zero-memory  nonlinearity  is  handled  as  a transforma- 
tion of  random  variables  on  the  incoming  signal  and  noise.  As  a 
result  of  these  nonlinearities,  the  output  envelope  pdf's  are 
finite,  i.e.  the  envelope  cannot  exceed  some  point  t^ . After 

more  transformation  of  variables,  we  arrive  at  pdf's  of  each  of 

t h 

the  y^ 's.  These  pdf's  are  also  finite.  Hence,  the  k moment  of 
y , is  of  the  form 

k k k 

E(yJ}  = / y^  p(ya)  dy£  (84) 

None  of  the  noise  models  known  to  the  author  to  date,  including 
those  in  this  study,  are  sufficiently  rich  in  analytical  useful- 
ness that  closed  form  solutions  exist  for  all  the  computations 
of  interest.  It  is  interesting  to  note  that  even  in  the  case 
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of  Gaussian  noise,  if  a nonlinearity  is  present,  integrals  such 
as  Eq.  84  cannot  always  be  solved  in  closed  form. 


5.2  Simulation  Error-Rate  Results 


In  this  section,  we  present  the  simulation  estimates 
of  the  probability  of  error  for  various  CPMFSK  modulations. 
Estimates  are  presented  for  CPMFSK  of  different  modulation  index 
M,  different  deviation  ratios  h,  different  ratios,  different 
receiver  types,  and  at  a range  of  signal- to-noise  ratios. 

Because  of  the  vast  number  of  variables  of  interest,  it  is  of 
course  implausible  to  display  or  evaluate  results  which  are 
exhaustive  over  the  range  of  these  variables.  Hence,  the 
emphasis  has  been  to  understand  trends  and  the  general  phenomena 
of  using  different  receiver  structures  to  process  CPMFSK  signals 
which  are  corrupted  by  impulsive  noise. 


From  Eqs . 5,  6,  and  10,  the  baseband  equivalent  signal 
plus  noise  is  during  the  i*"*1  symbol  interval 


r ( t) 


d . Trh  ( t-  ( i — 1 ) T i-1 

exp  {j  [— = — + Tih  Z 

s j = l 


j<Mt) 

+ n (t)  e 


No  filtering  effects  have  been  included  in  the  results  to  be 
presented  in  this  section,  however  the  received  signal  plus  noise 
is  sampled  four  times  per  bit  (for  quaternary  modulation,  8 times 
per  symbol) . This  implies  an  RF/IF  bandwidth  (one-sided)  of  2/T^ 
Hz,  where  is  the  bit  period.  Reference  to  section  3.1 
indicates  that  filters  with  bandwidths  of  this  magnitude  will 
have  negligible  distortion  effects  on  the  signal,  hence  the 
filtering  effects  were  not  included  in  the  simulation. 
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During  each  iteration  within  the  simulation,  n random 

data  symbols  are  generated  which  correspond  to  d^,  . . . , dn . A 

complex  signal  array  is  then  constructed  corresponding  to  the 

first  term  on  the  r.h.s.  of  Eq . 85.  The  length  of  the  array  is 

4n  for  binary  modulation  and  8n  for  quaternary  modulation  (4 

samples  per  bit  X log2  M x n bits).  To  this  array  is  added 

complex  noise  samples  corresponding  to  the  second  term  in  Eq.  85. 

These  noise  samples  are  derived  from  the  simulation  noise  model 

2 

described  in  section  4.2  where  the  variance  of  the  noise  is  o 
2 

(Eg/ o ) l s an  input  parameter  to  the  simulation).  The  phase 
samples  of  the  noise  (f>  ( t ^ ) are  uniformly  distributed  on  the 
range  ( 0 , 2 ti  ) . 

After  the  signal  plus  noise  is  generated,  it  is  either 
processed  in  a linear  receiver  as  shown  in  Figure  1,  or  it  is 
first  passed  through  one  of  the  nonlinearities  of  Figure  3 and 
then  processed  in  the  linear  receiver.  The  linear  receiver  is 
also  a sampled  version  so  that  the  various  voltage  outputs  can 
be  modeled  as  a sum  rather  than  the  integral  of  Eq . 19.  That  is 

4 log2Mn 

x 'Re  E z ( t,  ) s * ( t,  , d , , D • ) (86) 

Aj  k=1  K K A ] 

The  Mn  different  outputs,  x^,  are  then  compared  to  make  a 
decision  on  d^.  The  number  of  errors  are  counted  and  the  simula- 
tion program  was  exercised  until  enough  iterations  were  taken  to 
assure  the  estimated  error  rate  was  reasonably  accurate.  To 
determine  this  number  of  iterations,  the  following  calculations 
were  used: 

Assume  errors  occur  with  a binomial  probability  p.  Hence, 
in  N trials,  the  expected  number  of  errors  is  Np  and  the  variance 
of  the  errors  would  be  Npq  where  q=l-p. 
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For  p very  small,  the  ratio  of  standard  deviation  to  mean  is 


/ Npq 
r Np 


/~Np 

Np 


(87) 


For  the  bulk  of  the  results  in  this  study,  it  was  desired  to  have 
an  error  bracket  e = i,  i.e.  the  simulated  points  are  accurate 
with  a one  standard  deviation  range  of  ± of  the  value.  From 
Eq.  87,  this  implies  that  Np  = 25.  For  an  estimated  error  rate 
of  p = 10-^,  up  to  25000  iterations  would  then  be  required.  This 
indicates  the  relatively  high  expense  of  determining  low  error 
rates  using  Monte  Carlo  simulation. 

Since  results  can  logically  be  split  up  depending  upon 
what  receiver  type  was  used,  the  following  sections  will  display 
results  for  linear,  clipper,  and  hole  punch  receivers  respectively. 


5.2.1  Simulation  Results  for  Linear  Receivers 

One  of  the  first  considerations  which  must  be  examined 
with  respect  to  CPMFSK  signaling  is  the  choice  of  the  deviation 
ratio,  h.  So  far,  no  general  techniques  have  arisen  for  determin- 
ing the  optimum  h,  and  different  approaches  have  been  used. 

The  value  of  h which  should  be  chosen  is  obviously  the 
one  that  minimizes  the  probability  of  error.  Some  authors  in  the 
past  have  alternately  examined  the  minimum  distance  criterion, 
wherein  a value  of  h is  chosen  which  maximizes  the  minimum  dis- 
tance between  pairs  of  signals  in  the  signal  set  [3,  24].  In 
this  study,  an  attempt  was  made  to  find  the  value  of  h for  which 
PrM(£)  was  smallest. 

Table  6 shows  the  error  rate  results  for  various  binary 
CPFSK  signals,  including  MSK  (h=l/2) . The  results  were  generated 
at  3 different  ratios,  however  at  each  ratio  only  one 
signal-to-noise  ratio  (E_/o^)  was  used. 

■D 
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Table  6:  Error  Rates  for  Binary  CPFSK  vs. 
Deviation  Ratios 


Vd  = 2 dB 

V = 8 dB 
a 

Vd  = 14  dB 

h 

En/o2  = 3 dB 

D 

Eg/a2  = 3 dB 
-2 

E /o2  = 3 dB 

-2 

.4 

2.45  x 10~2 

-2 

2.13  x 10 

-2 

1.26  x 10 

-2 

. 5 

2.44  x 10 

-2 

1.96  x 10 

„ -2 

1.30  x 10 

-2 

2/3 

1.5  x 10 

1.53  x 10 

1.1  x 10 

.715 

1.46  x 10“2 

1.36  x 10-2 

1.06  x 10-2 

. 75 

1.48  x 10~2 

1.46  x 10"2 

1.1  x 10-2 

.8 

1.77  x 10~2 

1. 50  x 10~2 

1.1  x 10"2 

All  the  data  except  for  h = 0.5  was  derived  with  a reciever  of 
memory  length  n = 3,  i.e.  three  bit  are  observed  before  a decision 
is  made. 


Three  things  must  be  noted  with  respect  to  the  data  of 

Table  6: 

(1)  The  data  is  derived  from  simulation  results,  so 
there  is  a certain  amount  of  expected  randomness 
in  it . 

(2)  The  error  rates  for  values  of  h other  than  those 
listed  is  not  necessarily  a smooth  function  which 
could  be  surmised  from  interpolation.  [3]  This  is 
why  the  data  was  presented  tabularly  rather  than 
in  a figure.  The  h's  which  were  investigated 
were  found  to  be  locally  optimum  points  on  the 


61 


Gaussian  channel,  but  this  statement  cannot 
be  made  in  general  on  the  impulsive  channel. 

2 

(3)  The  signal-to-noise  ratio  (E D/a  ) was  kept 

D 

constant,  so  it  is  conceivable  that  other 
values  of  h could  be  better  at  different 
Eb/o2 

Nonetheless,  the  data  does  indicate  that  h = .715  (which  was 
found  optimum  in  the  Gaussian  case)  is  a good  choice  for  binary 
communications  when  error  rate  is  the  primary  consideration. 

Table  7 shows  similar  data  for  quaternary  CPMFSK  signaling. 
Again  all  data  was  derived  using  n = 3 which  implies  a receiver 
with  memory  extending  over  3 symbols. 

Table  7:  Error  Rates  for  Quaternary  CPMFSK  vs.  Deviation  Ratio 


V,  = 2 dB 
d 

V-.  = 8 dB 
d 

Vd  = 14  dB 

h 

E / Q2  = 0 dB 

2 

En/a  = 0 dB 

E /a2  = -3  dB 

.4 

6.3  x 10  2 

3.33  x 10-2 

3.47  x 10"2 

-2 

-2 

-2 

. 5 

6.4  x 10 

3.87  x 10 

3.87  x 10 

.6 

4.8  x 10"2 

3.73  x 10-2 

4.0  x 10“2 

.8 

4.53  x 10"2 

3.06  x 10'2 

2.94  x 10-2 

.9 

6.67  x 10“2 

3.33  x 10-2 

3.07  x 10-2 

The  same  three  cautions  given  above  also  of  course  apply  to 
Table  7.  h values  of  .4,  .6,  and  .8  were  shown  to  achieve  good 
performance  improvement  on  the  Gaussian  channel.  h = .5  is 
minimal  orthogonal  spacing.  From  Table  7,  it  appears  that  h = .4 
and  .8  still  achieve  error  rates  lower  than  orthogonal  signaling, 
however  h = .6  does  not  seem  to  be  as  good  for  Vd  = 14.  This 
could  be  a statistical  anomaly. 
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It  was  determined  that  for  linear  receivers,  the  best 
signaling  scheme  in  terms  of  probability  of  error  is  quaternary 
CPMFSK  with  h = .8.  This  agrees  with  previous  Gaussian  results  [4] . 
All  of  the  other  binary  and  quaternary  signals  which  achieved 
better  error  rates  than  MSK  were  somewhere  between  MSK  and  quater- 
nary CPMFSK  with  h = .8.  Consequently,  all  the  data  on  linear 
receivers  and  most  of  the  data  on  nonlinear  receivers  will  be 
presented  in  terms  of  MSK  performance  and  quaternary  CPMFSK  with 
h = .8  performance.  The  MSK  performance  can  be  considered  a 
benchmark  because  it  has  been  used  on  some  impulsive  channels 
( VLF/LF ) and  it  is  the  minimum  orthogonal  frequency  spacing.  The 

quaternary  CPMFSK  with  h = .8  signaling  is  the  best  that  has  been 

2 

found  so  far  in  terms  of  error  rate  vs.  E Q/o  . It  is  probably  a 
reasonable  assumption  to  imagine  that  all  the  other  modulations 
which  are  better  than  MSK  would  achieve  error  rates  between  MSK 
and  quaternary  CPMFSK  with  h = .8. 

A curve  showing  the  linear  receiver  error  rate  per- 
formance of  MSK  and  quaternary  CPMFSK  with  h = .8  for  three 
different  ratios  is  shown  in  Figure  20.  The  degradation 
caused  by  impulsive  noise  is  apparent  from  the  flattenning  of 
the  curves  as  the  noise  becomes  more  impulsive  (higher  ratio) . 

In  Figure  20  the  probability  of  bit  error  is  plotted  with 
2 

respect  to  Eg/a  where  Eg  is  the  received  energy  per  bit  and 

is  the  variance  of  the  impuslive  noise.  For  M-ary  modulation, 

the  received  energy  per  symbol  E is  related  to  ED  by 

S B 


Es  = 1o^2  M'EB 


(88) 


Again,  it  is  reiterated  that  Figure  20  is  the  result  of  Monte 

Carlo  simulation,  so  that  plotted  points  are  expected  to  have 

some  statistical  jitter  on  them.  Nonetheless,  it  is  seen  that 

2 

CPMFSK  can  offer  improvement  in  En/o  from  1.5  dB  to  3.5  dB  at 

_ 3 B 

error  rates  of  10  . This  corresponds  to  the  2.3  dB  improvement 
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Figure  20:  Linear  Receiver  Performance 


6 


possible  on  the  Gaussian  channel  [4]  at  the  same  error  rate.  On 
the  other  hand,  as  mentioned  above,  the  overall  performance  of 
the  linear  receiver  is  much  worse  when  impulsive  noise  rather 
then  Gaussian  noise  is  present.  This  was  to  be  expected  from 
theoretical  studies  [12,  16,  25,  26]  which  derived  optimum  re- 
ceiver structures.  Although  the  details  of  these  structures 
depend  upon  knowledge  of  the  noise  pdf,  it  is  clear  that  for 
impulsive  noise,  linear  receivers  are  decidedly  suboptimum. 

5.2.2  Simulation  Results  for  Nonlinear  Receivers 

In  addition  to  the  considerations  given  in  the  last 
section  on  linear  receivers,  additional  related  considerations 
must  be  addressed  in  the  case  of  nonlinear  receivers.  Only 
nonlinear  receivers  which  are  composed  of  the  linear  receiver 
of  Figure  1 preceded  by  zero-memory  nonlinearities  are  considered 
in  this  study.  Furthermore,  only  two  types  of  nonlinearities 
are  considered,  namely  the  hole  punch  (or  blanker)  and  the 
clipper.  The  transfer  functions  of  these  are  shown  in  Figure  3. 
For  both  of  these  nonlinearities,  one  further  parameter  to  in- 
vestigate is  the  threshold  level  t^.  All  results  in  this  study 
were  determined  where  the  threshold  level  is  set  relative  to 
the  mean  of  the  signal  plus  noise  envelope.  As  indicated  in 
Eq.  9,  the  received  signal  preceding  the  nonlinearity  can  be 
represented  as 

r ( t)  = Re  | v(t)e^a(t)e:’U)ct|  (89) 

The  threshold  level  can  then  be  represented  as 


OO 

t^  = k / vp(v)  dv  = kv  (90) 

a 

where  p(v)is  the  pdf  of  the  signal  plus  noise  envelope,  and  k 
is  a proprotionality  constant.  Results  are  displayed  with  t^ 
indicated  in  dB  relative  to  v. 
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Figure  21  shows  a plot  of  v for  three  different  ratios 
and  a normalized  signal  energy  Eg  = 1 . The  simulation  noise  model 

explained  in  Section  4.2  was  used  and  results  are  displayed  as  a 

2 2 - 
function  of  Eg /a  . Clearly,  as  Eg/a  gets  larger,  v approaches  1 

as  expected.  The  rest  of  this  section  presents  results  for  the 

clipper  and  hole  punch  nonlinear  receivers. 

5. 2. 2.1  Clipper  Results 

The  nonlinear  clipper  has  an  envelope  transfer  function 
as  shown  in  Figure  3a.  It  is  perhaps  the  most  commonly  used  non- 
linearity for  noise  suppression.  The  technique  shown  in  Figures 

22  and  23  was  used  to  determine  the  optimum  threshold.  The  sig- 

2 

naling  waveform,  noise  V,  ratio,  and  E /a  were  all  fixed  while 

Q D 

the  nonlinearity  threshold  t^  was  varied.  Threshold  values  of  0, 

+1,  +2,  +3,  +6,  and  +10  dB  relative  to  v were  investigated.  Al- 
2 

though  E_/a  remained  fixed  for  each  curve,  as  can  be  seen  in  the 
® 2 

figure,  different  Eg/a  were  used  for  different  signals  and 
ratios.  This  was  done  to  more  clearly  display  the  effect  of 
varying  t^.  As  is  seen  in  Figure  22,  the  best  threshold  setting 
in  the  investigated  range  for  the  MSK  signals  is  -10  dB  relative 
to  v.  On  the  other  hand,  in  Figure  23,  the  curves  for  = 2 dB 
and  = 8 dB  show  a clear  improvement  for  t^  = -2  dB  and  -6  dB 
relative  to  v respectively.  It  is  suspected  that  optimum  threshold 
values  might  exit  below  -10  dB  for  the  other  signals  and  ratio, 
but  these  were  not  investigated.  For  each  of  the  points  in 
Figures  22  and  23,  the  Monte  Carlo  noise  model  was  exercised  from 
the  same  starting  number.  Hence,  the  relative  differences  be- 
tween points  should  be  fairly  accurate. 

Figure  24  shows  the  probability  of  bit  error  when  the 
optimum  clipper  threshold  was  chosen.  The  word  'optimum'  is  used 
where  it  is  understood  that  the  threshold  setting  is  in  the 
range  as  shown  in  Figures  22  and  23.  A number  of  things  become 
clear  when  investigating  Figure  24.  For  one  thing,  comparison 
with  Figure  20  shows  that  significant  performance  improvement 
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Figure  23:  Clipper  Performance  vs.  Threshold 
Level  (Quarternary) 


Figure  24:  Optimum  Clipper  Performance 
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over  linear  receivers  can  be  achieved  when  nonlinear  noise 

2 

suppression  techniques  are  used.  The  improvement  in  E_,/o 

13 

is  on  the  order  of  1.5  dB  when  the  V.  ratio  is  2 dB , however 

Q 

it  is  on  the  order  of  21  -22  dB  when  the  V.  ratio  is  14  dB . 

d-i 

These  improvement  are  for  error  rates  of  10  . 

It  is  also  clear  from  Figure  24  that  CPMFSK  still  achieves 

2 

a lower  error  rate  thar.  MSK  at  the  same  Eg/o  , however  the  per- 
formance improvement  does  not  seem  to  be  as  great  when  nonlinear 

receivers  are  used  as  when  linear  receivers  are  used.  Not  only 
2 

is  the  Eg/a  improvement  of  CPMFSK  over  MSK  not  as  great  in 
Figure  24  as  it  is  in  Figure  20,  but  the  improvement  seems  to 

become  less  as  the  noise  becomes  more  impulsive.  For  example, 

- 3 2 

at  error  rates  of  10  , CPMFSK  offers  an  advantage  in  E /o  over 

B 

MSK  of  dB  for  = 2 dB,  whereas  the  improvement  appears  to 

be  about  0.5  dB  for  = 14  dB . On  the  other  hand,  it  must  be 

remembered  that  all  the  points  in  Figure  24  are  derived  from 

simulation  so  there  is  a certain  expected  randomness  in  their 

accuracy.  Also,  simulation  expenses  precluded  the  determination 

of  points  much  below  10  ^ , so  it  is  unclear  if  the  V = 14  dB 

d 

curves  will  tend  to  diverge  at  low  error  rates  as  it  appears  the 
lower  Vd  curves  do. 

The  importance  of  picking  the  correct  threshold  level  t^ 
is  again  shown  in  Figure  25,  where  all  results  are  shown  for  a 
threshold  setting  t^  = v (0  dB) . Comparison  with  Figure  24 
shows  that  for  ratios  of  8 and  14  dB,  the  optimum  threshold 

setting  achieves  perforance  improvements  (required  ER/a  for 

- 3 b 

bit  error  rates  of  10  ) of  2-3  dB  over  what  would  be  achieved 

if  t^  were  simply  chosen  to  be  v. 

5. 2. 2. 2 Hole  Punch  Results 

The  nonlinear  hole  puncher  has  an  evelope  transfer 
function  as  shown  in  Figure  3b.  As  in  the  case  of  the  clipper, 
an  investigation  of  the  optimum  threshold  setting  t^  was  done  for 
the  hole  punch.  Results  are  shown  in  Figures  26  and  27.  In 
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Figure  27:  Hole  Punch  Performance  vs.  Threshold  Setting 
(Quaternary) 


contradistinction  to  the  variation  of  clipper  error  rate  vs.  t^, 
Figures  26  and  27  show  that  hole  punch  performance  is  very 
sensitive  to  t^  and  optimum  values  of  t^  are  more  evident. 

Also,  all  of  the  optimum  values  of  t^  are  greater  than  or 
equal  to  -1  dB  relative  to  v.  This  is  intuitively  satisfying 
because  it  indicates  that  performance  degrades  if  the  threshold 
value  is  set  at  a value  lower  than  the  desired  signal  level. 

Figure  28  shows  the  error  rate  results  when  the  optimum 

threshold  setting  is  used.  Again,  it  is  clear  that  nonlinear 

processing  of  the  signals  can  achieve  significantly  improved 

performance  over  the  linear  receiver  (refer  to  Figure  20) . 

On  the  other  hand,  comparison  of  Figures  24  and  28  shows  that 

optimum  hole  punch  performance  is  not  as  good  as  optimum  clipper 

2 

performance,  where  again  the  E /a  required  to  achieve  an 
- 3 “ 

error  rate  of  10  is  used  as  a criterion  for  goodness. 

Figure  29  shows  the  error  rate  curve  if  the  threshold 
value  t^  is  set  equal  to  v.  Since  the  optimum  thresholds  are 
close  to  v from  Figures  26  and  27  for  all  curves  except  = 2 dB, 
the  error  rate  curves  are  not  significantly  different  from 
Figure  28.  Comparison  with  Figure  25  shows  that  the  hole  puncher 
with  t^  = v outperforms  the  clipper  with  t^  = v for  = 8 
and  14  dB  impulsive  noise  ratios.  This  is  pontentially  useful 
in  a situation  where  an  AGC-type  signal  can  estimate  v but  little 
or  no  knowledge  is  available  about  the  Vd  ratio.  In  such  a 
situation,  it  might  be  desirable  to  set  the  threshold  level 
approximately  equal  to  v and  a hole  punch  can  give  better  per- 
formance than  a clipper  for  largely  impulsive  noise  (V^  ratio 
high)  . 

Both  Figures  28  and  29  show  that  CPMFSK  can  still  out- 
perform MSK  when  hole  punchers , are  used.  However,  there  is  a 
tendency  for  performance  improvement  of  CPMFSK  over  MSK  to  become 
less  pronounced  as  the  noise  becomes  more  impulsive  (higher  V^) . 
This  was  also  seen  in  the  case  of  noise  suppression  by  clipping 
(Section  5 . 2 . 2 . 1 ) . 


• MSK,  Vd  =2  dB,  THRESHOLD  = +6  dB 

▲ MSK,  V,  8 dB,  THRESHOLD  =+l  dB 

d 

■ MSK,  Vj  = 14dB,  THRESHOLD  = 0 dB 

d 

O QUATERNARY  CPFSK,  n=3,  h=.8,  Vd  = 2 dB,  THRESHOLD  = +6  dB 

A QUATERNARY  CPFSK,  n=3,  h=.8,  = 8 dB,  THRESHOLD  = + 1 dB 

Q QUATERNARY  CPFSK,  n=3,  h=.8,  = 14  dB,  THRESHOLD  = -)  dB 
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Figure  28:  Optimum  Hole  Punch  Performance 
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5.3  Summary  of  Error  Rate  Results 

At  the  present  time  quantitative  error  rate  results  exist 
only  using  Monte  Carlo  simulation,  however  in  all  cases  which 

2 

were  investigated,  a general  CPMFSK  signal  required  less  E n/o 

_ 3 " 

to  achieve  an  error  rate  of  10  than  MSK . Significant  per- 
formance improvement  is  possible  if  nonlinear  noise  suppression 
devices  are  used  preceding  a linear  receiver.  If  optimum 
thresholds  are  used,  bandpass  clippers  achieve  better  performance 
than  bandpass  hole  punches,  however  if  thresholds  are  chosen  to 
be  equal  to  the  mean  value  of  signal  plus  noise  envelope,  the 
hole  puncher  achieves  better  performance. 
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6.0  Power  Spectra  for  CPMFSK 

The  one-sided  power  spectral  density  for  CPMFSK  is  known 
from  [1]  and  can  be  written  in  a normalized  form  as 
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where 


A = 


6 = 


(“  " wc)Ts 

2i  1 


d h 

Yn  = (6  " • 

B = Cos(Yn  + Ym}  ~ C Cos(Yn  + Ym  ~ 27Tg)  , 
1 + C2  - 2C  Cos  2tt 6 

and 

2 M/2 

C = ^ £ CosTTh  ( 2n  - 1) 

n=l 


(92) 


This  is  a formidable  looking  set  of  expressions,  but  they  are 
relatively  easy  to  calculate  on  a computer.  Figure  30  shows  the 
results  for  the  special  case  of  MSK  (h  = .5)  whereas  Figures  31 
and  32  show  the  results  for  binary  and  quaternary  modulation. 
Also  indicated  on  the  figures  are  the  numerically  determined 
50  percent  and  99  percent  normalized  bandwidth. 


It  must  be  noted  that  Figure  32  is  plotted  in  terms  of 

the  symbol  time  T and  symbol  energy  E whereas  Figures  30  and 

s s 

31  are  plotted  in  terms  of  the  bit  time  TD  and  bit  energy  En . 

D D 
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Figure  30:  One-Sided  Spectrum  of  MSK 


Normalized  Binary  CPFSK  One-Sided  Power  Spectral  Density 


50%  BW  99%  BW 


Figure  32:  Normalized  Quaternary  CPFSK  One-Sided  Power  Spectral  Density 


The  two  are  related  by  Eq.  88  and 


Tg  = log 2 M-Tb  (93) 

so  that  if  Figure  32  were  plotted  in  terms  of  T_,  and  ED,  the 
curves  would  appear  twice  as  high  and  half  as  narrow  as  they 
do  in  Figure  32.  The  quoted  50  percent  and  99  percent  band- 
widths  of  Figure  32  would  also  be  halved.  This  is  shown  in 
Figure  33  which  displays  again  the  spectra  of  MSK  and  quaternary 
CPMFSK  with  h = .4,  where  both  curves  are  plotted  on  the  same 
scale.  It  is  interesting  to  note  that  the  quaternary  CPMFSK 
with  h = 0.4  has  a smaller  99  percent  bandwidth  than  MSK. 
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Figure  33:  Normalized  Spectra  of  MSK  and  Quaternary  CPMFSK  (h 


7. 


Conclusions  and  Recommendations  for  Further  Work 


The  results  of  this  study  show  conclusively  that  the  use  of 
CPMFSK  as  a modulation  on  channels  with  additive  atmospheric 
noise  requires  less  signal  to  noise  ratio  to  achieve  a given 
error  rate  than  would  be  necessary  with  more  conventional  modu- 
lations such  as  PSK,  MSK , or  CFSK.  Furthermore,  transmitted 
signal  power  spectral  densities  have  comparable  or  smaller  99% 
bandwidths  than  these  more  conventional  modulations. 

A Monte  Carlo  simulation  to  determine  error  rate  performance 
was  developed  and  exercised.  A theoretical  approach  to  determine 
error  rate  based  upon  a Gram-Charlier  series  was  formulated, 
however  time  constraints  precluded  its  full  development  and  ap- 
praisal . 

Specific  recommendations  for  further  work  include  the  fol- 
lowing : 

1.  The  theoretical  approach  should  be  further  investigated 
to  determine  its  applicability  and  feasibility.  If  the  approach 
proves  to  be  feasible,  the  following  subjects  should  be  investi- 
gated . 

a)  Determine  the  optimum  deviation  ratios  to  ascer- 
tain if  they  differ  from  those  which  are  optimum  in  Gaussian 
noise. 

b)  Determine  the  effect  of  correlated  noise  on  the 
error  rate  performance. 

2.  Other  additive  atmospheric  noise  models  may  be  more 
amenable  to  either  analytic  determination  of  required  parameters 
or  to  quicker  numerical  determinations.  This  should  be  investi- 
gated . 

3.  In  the  HF  band,  fading  becomes  an  important  signal  im- 
pairment along  with  additive  atmospheric  noise.  These  fading  ef- 
fects should  be  included  in  both  the  theoretical  and  simulation 
phases  of  the  study. 
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4.  A hardware  development  to  determine  implementation  and 
quantization  losses  should  be  instigated. 
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